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I. 


Ein großes Interesse liegt, wie bekannt, im Studium der Bewegung eines elektrisierten Teil- 
chens in einem elektrischen und in einem demselben überlagerten magnetischen Felde oder in 
einem elektromagnetischen Feldet). In dieser Mitteilung werde ich über die Resultate einiger neuer 
Untersuchungen berichten, die die Bewegung eines Elektrons in einem stationären magnetischen 
Felde H und in einem demselben überlagerten stationären elektrischen Felde E betreffen. Das 
magnetische Feld ist (im ersten betrachteten Fall) zu der x-Achse symmetrisch und von einem 
Potential @(x, y) abhängig; das elektrische Feld ist nur als Funktion von x oder als ein konstanter 
Vektor gedacht. Bezeichnen wir mit © den von der Bewegungsebene des Elektrons und von einer 
fixierten Ebene gebildeten Dyeder (mit der Kante x). Die relativistische Differentialgleichung 


der Elektronenbewegung ist: 
1 


md) =eE+ind, m-mda—e ?|, 


wo d die Geschwindigkeit des Teilchens, e die Ladung des Elektrons, m, die Ruhmasse und c die 
kritische Geschwindigkeit sind. Außerdem haben wir: 


re plz, Y) Ye, y) _ 9 Y) 
dx rg Eee dy ro’ 


wo y(z, y) die „Funktion von Beltrami‘ und g ein konstanter Parameter sind. Die Bewegungs- 
differentialgleichung verwandelt sich in folgender Weise: 


H = grad pa, 3), 


Ge ge — » 5 

md) = eE@ I +ulögrady—ygradQ) :........ (1), 
— ed al 
[/ = Einheitsvektor in der x-Richtung, © = 2 HAN Fr’ u= r = const]. 


Für die Differentialgleichung (1) gilt das erste Integral der Energie in relativistischer Form: 


2 2 2232| 2 
TER dr = eu = me | TU nl = h, (K=cons) 2.2.1202), 


c2 


ech x 
Multiplizieren wir (1) skalar mit nn so haben wir: 


00 


P 
= (m y? 6) = u (O grad y — ygrad ©) io=- ud, d.h. 


mPÖ+uy=C=cons . EEE EBENE 72 (0) 
ein neues erstes Integral?). Jetzt beachten wir, daß 


dere op "id N a 2 2 
zu) md: mv) x 


BT me , 
du = my)—my® 


*) Vortrag, gehalten auf der wissenschaftlichen Jahrestagung derGAMM in Hannover, 19.—23. Mai 1959. 


1) A. Pignedoli, Zeitschrift für Angewandte Mathematik und Mechanik 88 (1958), S. 343—348. 


2) Für das Studium der Bewegung eines Elektrons in einem zu einer Achse symmetrischen magnetischen 
Felde in der klassischen Mechanik vgl. C. Agostinelli, Atti della R. Accademia delle Scienze di Torino, vol. 74, 


1938—39. 
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ist. So haben wir, für die betrachtete Bewegung, das System von Differentialgleichungen: 


en uy.d, 2 (CERWiT ER Tun 
FRUDERP ORT m y? ’ my) = my? T 0y my? Zn: (4). 


Führen wir in die Bewegungsdifferentialgleichung den Zeitparameter 7, der mit £ durch 
1 


dt _ m u 1 -=1—-@+P+% le] * a (5) 


de m y1— ve 


verbunden ist, so wird die Bewegungsdifferentialgleichung die folgende: 


ar = dP_ + 
ma memEtem , NH 2 gn u a Na a (6), 
und für diese gilt das erste Integral: 
m=m+eU()/ce, (m, = const)" „uw. 2 Se er (A). 
Infolgedessen können wir das System (4) schreiben: 
d 0 i 
mm). E + | 
(mM; 
ET E le 1 (C— uy) 
er, TE Ni 
daraus folgt das erste Integral: 
i 3 C— 2 
m: (+1) =2em, U) + = U?(x) er +28, (Kenn Fe 
Ist das elektrische Feld ein konstanter Vektor E=E- vB so haben wir: 
2 2 
ne +) = Lem Er+°, He (C er 4 2K, 
und, wenn E = 0 ist, dann ist die Masse m konstant (m = m,), und es gilt das erste Integral: 
er 2 « 
nee Sr a 7 


y° 
R Hat das magnetische Feld eine zur x-Achse orthogonale Symmetrie-Ebene, so haben wir 
&4 =( und folglich: 


Ge 
Ys 
dı 
t—I,=m | — = 2 er er 8 
i Az, - (C— u yp(0, yYP 
Y 
Yı 
. [CE — u y(0, y)] dı 
9— 09 eu ale a yl 1 dy be 
yy2K y? — IC — u y(O, y)? 
v 
Il. 
Be dene 3 a PLANE RER I nn 
un studieren wir den Fall, in welchem E, = E,(x) = E(x) = a Eh =-E=0, 
V 
H,=- 4,0) ed ®e #0, H=H,=D0 ist. 


Die Bewegungsdifferentialgleichungen lauten: 


d ; R d k f 
Fr (mx) =eE(x) +eyH(«), Fri (my) =—ezxH(), e {m 2),==0 Wezrldn 
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Man kann dieses System integrieren und findet: 


- Ua m) ar 
ee - | TEA 
V: UX2) + 2em, U) +2eq, Va) —eV%a) +2h, 


y—n=H [9, — e V(@)] dx 
Vs U?(x) SR 2e m, U(x) + 2e go V(x) —. e2 V:(x) + 2h, 


AED): 


Ze er we (13). 
z UX2) + 2em, Ula) +2eg Ve) —e Va) + 2h, 
(do Fo, A, const) 


Der besondere nicht-relativistische Fall, in welchem m = m, und 7 = t ist, interessiert die 
pinching-Bewegung eines Plasmas. In diesem Falle haben wir: 


7 — = | dx = a VE RE RE (14), 
m,& = eE(x) + ey, H(x) — H6) DE Ar ehe er (19) 


und infolgedessen: 


% 


De : = er (16), 
; 2 E r 
+ V: + =|| E(&) + ey, H(&) te ae] dx 
NN er en ————————— 
+ n || E(x) + ey, H(&) — 1, H@) ze] dx 
x et) I a ychrr 


mu|/ a3 + 2 | Ä El) + ei. H@)— © Ha) 20) da 
0 0 


%o 


III. 


Jetzt betrachten wir (vgl. Bild) dieBewegung eines relativistischen Elektrons P in dem magne- 
tischen Felde einer kreisförmigen Windung vom Radius a, in welcher ein ununterbrochener elek- 
trischer Strom I fließt, und in einem demselben überlagerten axialen konstanten elektrischen 
Felde E. 


Für die Bewegungsdifferentialgleichung des Elektrons 


md) =eE—e6grad Inlap [er ney nase 


v 


+egradO- 5, 2 to [er ney ne 
) 


(18), 
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wo man mit J, die Besselsche Funktion erster Art und erster Ordnung bezeichnet, haben wir 
das erste Integral der Energie 


m mt +:gz, (mf. == const). ; 2 KL ae ae 
und das neue erste Integral 
E oo 
m®O—2nlaeofer"J,l(osJlas)d=l=cont...... (20). 
0 "u 
Die Bewegungsdifferentialgleichung können wir also schreiben: 


|c+ Zara] BEAT, [OrAcer] 


m 0°? 


mn) eure grad (Oxtap ferera, (os) J, (ads) 


Hemd 0. ar tee e nes ER De (21), 
) 


und das Problem reduziert sich auf das Studium der Bewegungen eines relativistischen Teilchens 
mit Maupertuis-Masse m in der Ebene (9, z) in einem Kraftfelde mit den Komponenten 


{6} 2 
es C+2rnlaeo[ e*"J (os) Jı(as)ds 
0 


am, Dur 


und 


13 eE 1 (CH2alaeet ETUI HBETE 
ern Fr * 2 ” 
m 02 eE|m sl 2 


Für diese Bewegungen haben wir das erste Integral: 


[0] 


I 


j 1 A 1 C+2nlaeo[e*::J (os) J,(as) ds 
ea * le ER 0 . Bi 1 
,(@ + 22) = eEim a 5 . + const } (22). 


Ist E =+ 0, so können wir solche Bewegungen in der Richtung der z-Achse betrachten. Für 
diese Bewegungen findet man: 


_[  @*+eEzje) dz AR u 
ee cYm* + eBze m" ae m. [gr —) me | (23). 


Im Falle E= 0, gibt es die Möglichkeit von Bewegungen in der Ebene z = 0. Es gilt die 


Differentialgleichung: 
LH Re 2. 
En (m 0) = — 5m ee - ): (= 0)7 . 7, GT are 


; le®./3 1, ala [8 1 2 
wo (po) = 0 für o = 0; Y(o) = 2 6.7.28 )ür <a; o(0) = — ee 2) 


[4 


für 0 >a; G = hypergeometrische Funktion von Gauß. In diesem Falle haben wir außerdem 
c? (m*2 — mi) | 


m = m* (= const), v2 
m’*2 


= vo, (= const). 


a  . 
SE en; = 
m? 0° 


. (26). 


Es folgt: 


+ m*o a [C- 7% -ere)P | 


 m*2 g2 


! 
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Für die Bestimmung von o im Falle von stationären gleichförmigen kreisförmigen Bewe- 
gungen (0 = 0, 0 = R), die möglich sind, haben wir die zwei transzendenten Gleichungen: 


mnoe—C+eoo)=0, m’nwoe+C-erd)=-0 ......M 


und man findet © immer aus der Gleichung mt ® 6 — e oe) = C. 


Endlich betrachten wi den Fall von Bewegungen in der Nähe der z-Achse der Windung. 
Für diese Bewegungen haben wir das System von Differentialgleichungen: 


d LEO NG, nela:o 
m, md))=—; Ah 4 ar] 
d EN 1#941G melao R 
my(m2)=eEm ER "aram| DU (28); 
a an 
ne Tara © | 


und es gilt das erste Integral der Energie: 


I 2 fde\e az 
5 7% (®) +z)]- one 


2.27, ET URN re Ne 
| p) B F (a + 2) as (a? + 22) Ei 5) em? — m) (29): 


Vernachlässigen wir in dem Potential U die Größen in o? und betrachten C als eine infinite- 
simale Konstante, so haben wir Bewegungen mit 0 = o, (infinitesimale Konstante). Für diese 
Bewegungen findet man: 


age 


” ce N zel®C 
eE \: RETTEN r 
ee | lee 


Betrachtet man in U die Größen in 0%, und nimmt man C = 0, so findet man die Differen- 
tialgleichungen: 


Be, dz (30). 


(a? En zy2 


d?o ea U 0 
de "m(@+z2%° 


deze 9 eE Lone ]° 080° 
ee * | —_.m2 ee 
N re „fez|m er =] 3 | 


. (831). 


Für diese Gleichungen gilt das erste Integral der Energie: 


do\2 dz\? EN g 202 Palo? | 
me|(%) +2) Im +3) - m) Eure ee 182): 


Sei M? = const größer als das Maximum 2? e? I® a®/my des Koeffizienten von o in der ersten 
Gleichung des Systems (31) und vergleichen wir diese Gleichung mit der Gleichung der schwin- 


2 
genden Bewegungen (33) - + M?o=0. Wenn o,(T) und o,(rT) zwei Integrale — mit Be- 


rücksichtigung von (31) und (33) — sind und wenn diese Integrale denselben Anfangsbedingungen 
für = r, genügen, so finden wir, daß o,(T)/oz(T) rechts von 7 = r,, eine zunehmende Funktion 
von r (mit Anfangswert 1) ist?). Jetzt studiere ich einige Klassen von Bewegungen in einem zeit- 
lich veränderlichen magnetischen Felde B und in einem ihm überlagerten elektrischen Felde E, 
das mit dem ersten durch die Maxwellschen Gleichungen verbunden ist. Ein besonderer Fall 
ist der Fall der Bewegungen, die bei den Beschleunigern der Physik auftreten‘). 


Anschrift: Prof. Dr. A. Pignedoli, Universita di Bologna, Istituto Matematico ‚Salvatore Pincherle‘‘, Bologna 
(Italien) 


3) Vgl. auch: A. Pignedoli, Atti della Reale Accademia delle Scienze di Torino, vol. 92 (1957—58). 
4) Vgl. A. Pignedoli, Rend. della Accademia delle Sceienze dell’Istituto di Bologna, 1959. 
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Ein allgemeines Maximalisierungverfahren 
Von Paul Stefan Pütter*) 


Es wird ein numerisches Verfahren beschrieben, um das Maximum einer Funktion zu finden, wenn Pr 
als Nebenbedingungen sowohl nr als auch Ungleichungen vorliegen. Das Wesen des Aachen JR 
liegt darin, einen Punkt unter Berücksichtigung der Nebenbedingungen in Richtung des Gradienten p 
Funktion, deren Maximum gesucht wird, zu verschieben. Das Verfahren wurde bereits für den Rechen- 
automaten IBM 650 programmiert und hat sich bisher gut bewährt. : - 

A numerical method is described for finding the maximum of a given function under a set of supplementary 
conditions containing both equations and inequalities. The essence of the method comsists in displacing a 
point in the direction of the yradient of the given function in such a way that the supplementary conditions are 
satisfied. The method has already been used for programming on the automatic computer IBM 650 and has 
been found to work very satisfactorily. £ j 

On decrit un procede numerique afin de trouver le maximum d’une fonction si des Equations aussi bien 
que des late ae presentes en tant que conditions secondaires. La nature de ce procede est de deplacer 
un point — en considerant les conditions secondaires — dans la direction du gradient de la fonction dont 
le maximum est recherche. Ce procede a dejü &t& pr&pare pour le calculateur automatique IBM 650 et a 
bien repondu & l’attente jusqu’ici. 

OnucbIBaeTcH YUCHEHHBLIÜ METON 1A OTIPeNeJIleHUA MAKCHMyMa PyHKIMM B CJIyyae, KOTAa 
B KayecTBe NOIIONEIHUTEJIBHEIX YCJIOBHÜ UMEWTCH KAK yPaBHeHHA TaK H HepaBeHcTBa. CyIlHOCTB 
IIpmeMa COCTOHT B IIepenBI;KeHHN TOYKH B HaIlpaBJIeHuNn TpayıueHnta Toü PyHKIIUM, MAKCHMyM 
KOTOpoü Tpeöyerca HauTH, IPNHyeM TpeöyercH COÖJMONEHNE NOMOJIHUTEILHLIX YCIOBHÜ. ITOT 
METON ÖbIJI 3AIPOTPaMMHPOBAH JIA ÖBICTPoNelcTBylINel 3AIEeKTPOHHOÄ BBIYHCHHTEIBHOÄ 
ManımHbı Tuma MbBM 650 u Has yHOoBJIETBOPHTENIBHLIE pe3yJIbTaTkl. 


1. Einführung 


Mit der Entwicklung der modernen elektronischen Rechenautomaten ist es möglich ge- 
worden, mathematische Probleme zu lösen, deren numerische Behandlung früher wegen des 
außerordentlichen Rechenaufwandes unmöglich war. Eine der wesentlichen Fragen, die die 
Praxis nun immer wieder stellt, ist die Frage nach dem Maximum einer Funktion unter be- 
stimmten Nebenbedingungen. Für den Fall analytisch faßbarer Funktionen liefert die Analysis | 
mit dem Lagrangeschen Multiplikatorenverfahren schon einen brauchbaren Weg, der zur Lö- | 
sung führt. Schwierig wird die Aufgabe aber, wenn die Funktionen zu einer allgemeinen Be- 
handlung zu unhandlich werden. 


Theoretisch läßt sich das Lagrangesche Multiplikatorenverfahren immer dann anwenden, 
wenn die Nebenbedingungen in der Form von Gleichungen gegeben sind. Zusätzlich treten 
aber bei fast jeder praktischen Aufgabe noch Ungleichungen auf, denen die Variablen genügen 
müssen. So muß zum Beispiel die Menge einer Substanz immer positiv sein. Maximalisierungs- 
aufgaben für den Fall, daß zu den Nebenbedingungen auch Ungleichungen gehören, wurden 
bisher fast nur im linearen Fall bearbeitet, und es gibt für den linearen Fall auch schon eine 
umfangreiche Literatur. Das im folgenden beschriebene Verfahren zeigt nun eine Möglichkeit, 
auch im nichtlinearen Fall derartige Aufgaben zu behandeln. 


‚Gesucht werden soll das Maximum einer Funktion / unter den Nebenbedingungen g, = 0, 
wobei gleichzeitig die Variablen noch bestimmten Ungleichungen genügen müssen.. Die Aufgabe 
läßt sich formulieren: 


Mu... 2) = Max; nu nee a (1), 
TER) = UN Beh... MSN) Te 
h-nz20 0 v1. N 0 
u a ee 


Die Funktionen werden als differenzierbar vorausgesetzt. Die Funktionen 9, sollen von- 
einander unabhängig sein. Die Ungleichungen (3) haben eine recht einfache Form, und es soll 


am Schluß noch kurz gezeigt werden, wie sich auch andere Ungleichungen in dieses Verfahren 
einbauen lassen. 


Das Ziel der vorliegenden Arbeit war, ein möglichst allgemein anwendbares numerisches 
Verfahren zur Ermittlung eines Maximums zu suchen, und es war daher zunächst die Frage zu 
klären, was man überhaupt von einem solchen Verfahren erwarten kann. Über die Funktionen f 
und g sollten neben der Differenzierbarkeit keine weiteren Voraussetzungen gemacht werden. 


*) Babcockwerke, Oberhausen, Zentrale Rechenabteilung. 
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‚Berücksichtigen muß man aber, daß ein Rechenautomat in der endlichen zur Verfügung stehenden 
Zeit nur endlich viele Aussagen über eine Funktion auswerten kann. Dann ist er aber grund- 
sätzlich niemals in der Lage, unter den angegebenen allgemeinen Voraussetzungen einen Punkt 
zu ermitteln, von dem feststeht, daß sich seine Koordinaten um nicht mehr als eine vorgegebene 
Strecke von den Koordinaten des Maximalwertes in dem betrachteten Bereich unterscheiden. 
Man muß daher bei einem numerischen Verfahren, bei dem man über die benutzten Funktionen 
keine näheren Kenntnisse besitzt, immer in Kauf nehmen, daß man möglicherweise in die Irre 
geführt wird. 

In der Praxis ist nun in den meisten Fällen bekannt, wo ungefähr das Maximum liegen 
kann. Nach den vorstehenden Überlegungen erscheint es daher zweckmäßig, sich darauf zu 
beschränken, einen vorgegebenen Ortsvektor r,, der die Bedingungen (2) und (3) erfüllt, solange 
in Richtung zunehmender /(r) zu verschieben, bis eine weitere Verschiebung keinen nennens- 
werten Erfolg mehr zu bringen scheint. Die Glaubwürdigkeit eines so ermittelten „Maximums“ 
läßt sich noch dadurch prüfen, daß man den Ausgangsvektor z, über einen großen Bereich 
variiert. Der hier beschriebene Weg wird im folgenden beschritten, um ein allgemeines Maximali- 
sierungsverfahren zu entwickeln. 


2. Maximalisierung ohne Berücksichtigung der Ungleichungen 


Es sei zunächst angenommen, daß die Ungleichungen nicht berücksichtigt werden müssen, 
das heißt, daß in (3) 
h, =+%, 18 ee) 


ist. Vom Ausgangspunkt r, ausgehend soll dann eine Folge £;, ta, ts gefunden werden, für die gilt 


IK) < Mi) <a) <--- a a Be 3) 
El aa ee a 


ER DE a a u EIN: 


Zu diesem Zweck wird zunächst r, um einen Vektor y d, verändert in 


ERENTO 

Dabei ist dv, ein Einheitsvektor, der folgende Bedingungen erfüllt: 
DV 9,0) = 0; a BE A a ren 
DRUM IDEEN REEL ERRA OS): 


Die Bestimmung der Konstanten y wird weiter unten noch angegeben. Die Bedingungen 
(7) sichern, daß durch v, infinitesimal im Punkte x, die Nebenbedingungen (2) erfüllt werden. 
Dann liefert (8) unter der Voraussetzung der Nebenbedingungen den Vektor des stärksten An- 
stiegs, das heißt, den Gradienten unter Berücksichtigung der Nebenbedingungen. Der Vektor b, 


hat die Form 
Des Al 2 AN On) re ae ld 


wobei die Konstanten }„ aus dem Gleichungssystem zu berechnen sind, das sich durch Ein- 
setzen von (9) in (7) ergibt. Die Größe A dient der Normierung. Da die Funktionen g,(£) un- 
abhängige Nebenbedingungen sein sollen, kann für die Praxis der Fall, daß in speziellen Punkten 
die Vektoren \V g, linear abhängig sein können, als nur mit sehr geringer Wahrscheinlichkeit 
realisiert unberücksichtigt bleiben. 

Daß der in (9) angegebene Vektor d, wirklich die Bedingung (8) erfüllt, sieht man leicht, 
wenn man ihn mit einem beliebigen Einheitsvektor tw, vergleicht, der auch (7) erfüllt. Zerlegt 
man nämlich tv, orthogonal in iv; und iv, wobei tv} und i/ beide (7) genügen und wobei gleich- 
zeitig iv, parallel zu d, ist, so folgt aus 


WE 
wegen (7) und (9), daß auch V / senkrecht zu w/ ist. Dann ist aber 
nVjzWmV/= DV: 
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025 na 0 m an SS2En nal a a ein u = 


x . . . # Krz f E 
Der in (6) neu gefundene Punkt; wird im allgemeinen die Nebenbedingungen (2) verletz 
Es ist daher notwendig, durch ein iteratives Verfahren die Nebenbedingungen wieder zu e 
füllen. Man kann dafür etwa den linearen Ansatz 57 


RE A Se 1 en ee . (10) E: 

benutzen, wenn man gleichzeitig für die Funktionen g,(£) die ersten Näherungen na ® 

| NEE EVEN) - - "rer 
einsetzt. Aus g, = 0 folgt dann: Sg 
0= 9, + Z0m V In) VA) ee (12). 


Das Gleichungssystem (12) liefert die &,, und aus (10) erhält man dann x, als nächste Näherung. 
Es ist solange zu iterieren, bis die Funktionen g, mit hinreichender Genauigkeit erfüllt sind. 
Wird die Rechnung mit einem Rechenautomaten durchgeführt, so ist es daher zweckmäßig, 
von Anfang an für jede der Funktionen g, eine Konstante a, anzugeben, die die erforderliche 
Genauigkeit bestimmt. Es ist dann zu iterieren, bis für alle u 


no a (13). 
Es bleibt nun noch die Größe y in (6) festzulegen. In dieser Festlegung liegt natürlich ad 
eine gewisse Willkür, und man wird sich vor allem von zwei Gesichtspunkten leiten lassen. - 


Wählt man y groß, so führt man die Funktion f in großen Schritten ihrem Maximum zu. Man | 
nimmt aber gleichzeitig in Kauf, daß durch die oben beschriebenen Iterationen nur recht lang- 
sam — vielleicht sogar niemals — die Bedingungen (13) erfüllt werden. Macht man y klein, so | 
wird das Maximum in vielen kleinen Schritten erreicht. 


Es wird daher folgendes Verfahren vorgeschlagen. Vor Beginn der Rechnung wird für 
jede Variable x, eine maximale Schrittweite s, vorgegeben, die vernünftig erscheint. Der Rechen- 
automat wählt dann zunächst 


en ae a ee EEE 


v 


wobei v, die »-te Komponente des in (7) und (8) definierten Vektors vb, ist und berechnet sich 
mit diesem y das neue x, aus (6). Dieses g, wird dann durch Iterationen verbessert, um die Neben- 
bedingungen zu erfüllen, und die Anzahl der Iterationsschritte wird gezählt. Es soll nun weiter 
von Anfang an eine Konstante n vorgegeben sein, die angibt, wieviel Iterationsschritte die Ma- 
schine maximal ausführen soll. Ist diese Zahl erreicht und noch keine Konvergenz erfolgt, so 
halbiert die Maschine sämtliche maximalen Schrittweiten s,, was auch für alle Zukunft gelten 
soll und versucht, mit dieser halbierten Schrittweite die Nebenbedingungen (2) bzw. (13) zu 
erfüllen. Die Maschine halbiert die Schrittweiten solange, bis entweder Konvergenz bei den 
anschließenden n Iterationen erfolgt oder eine untere Schranke s für die Schrittweiten erreicht 
ist. Haben die Schrittweiten die untere Schranke s erreicht, so hält die Maschine an. Erfolgt 
Konvergenz bei den in (10), (11), (12) beschriebenen Iterationen, so wird der neu ermittelte 
Punkt als r, angesehen, und die Maschine prüft, ob 


[> Mt). u: Dann 2 


Ist die Ungleichung erfüllt, so gibt die Maschine dem Bediener einige interessierende Daten — 
wie r, und /(g,) — an und fährt mit r, als neuem Ausgangspunkt in der Rechnung fort, das heißt, 
sie berechnet sich zunächst einen neuen Vektor v,. Ist die Ungleichung (15) nicht erfüllt, so 
startet die Maschine vom Punkt x, mit halbierter Schrittweite. Das gesamte Maximalisierungs- 
verfahren wird durchgeführt, bis die Maschine ihre Schrittweiten bis zur unteren Schranke s 
verkleinert hat, da man dann in den Schrittweiten so etwas wie eine Genauigkeit der Koordinaten 
des Maximums sehen darf. 


Einige Worte sind noch zur Differentiation zu sagen. Man wird zweckmäßig von Anfang 
an für jede Variable x, eine Konstante k, angeben und dann unter der v-ten Komponente des 
Gradienten einer Funktion 9 die Größe 


1 
Vıa=-EWa.. 3 three] nd 


verstehen. Die k, können für die einzelnen Variablen je nach der Stärke der Änderung der 
Funktionen in bezug auf diese Variablen um Größenordnungen verschieden sein 
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3. Einfaches Verfahren zur Berücksichtigung der Ungleichungen 


Wesentlich erschwert wird die gesamte Maximalisierung, wenn Ungleichungen als Neben- 
bedingungen zu berücksichtigen sind. Da derartige Nebenbedingungen aber bei fast jeder 
Problemstellung der Praxis die Lage des Maximalwertes entscheidend beeinflussen, ist es un- 
erläßlich, sich mit ihnen eingehend zu befassen. Es soll in diesem Abschnitt zunächst ein nicht 
einwandfreies, dafür aber recht einfaches Verfahren zur Einführung von Nebenbedingungen der 
Form (3) besprochen werden. 


Die Auswirkung von Ungleichungen als Nebenbedingungen besteht darin, daß aus dem 
vorher unendlichen N-dimensionalen Raum der Variablen ein in gewissen — oder allen — Rich- 
tungen beschränkter Teil als allein zulässig herausgetrennt ist. Dieser herausgetrennte Teil ist 
aber noch N-dimensional, da eine Verminderung der Dimensionszahl nur durch Ah, = j, erfolgt 
und daher als trivialer Fall ausgeschlossen werden kann. Das im vorhergehenden Abschnitt 
beschriebene Verfahren zur Maximalisierung kann daher ohne Beeinträchtigung solange fort- 
geführt werden, bis man erstmalig auf eine „‚Grenzfläche‘“ trifft, das heißt, bis eine Ungleichung 
verletzt wird. - 

Es sei nun angenommen, daß ein Ausgangspunkt r, vorliegt, der vollkommen im erlaubten 
Bereich liegt. Er soll auch auf keiner der Grenzflächen liegen. Verletzt dann der Punkt x, in 
(6) eine oder mehrere der Ungleichungen (3), so wird y in ein y* so verkleinert, daß x; keine 
Ungleichung verletzt und auf wenigstens einer Grenzfläche liegt. Bei dem anschließenden 
Iterationsverfahren zur Erfüllung von g, = 0 wird der Punkt auf diesen Grenzflächen festge- 
halten. Dadurch wird die Anzahl der Nebenbedingungen um die Anzahl der Grenzflächen, auf 
die der Punkt gekommen ist, erhöht. Die Grenzflächen stellen dann auch Bedingungen dar, die 
von den g, unabhängig sind. Denn wegen (7) ist d, senkrecht’zu allen V 9,, und andererseits 
hat v, eine Komponente in Richtung der Normalen der Grenzfläche. Dann ist aber die jetzt 
als verbotene Richtung hinzugekommene Normale der Grenzfläche linear unabhängig von den 
bisherigen verbotenen Richtungen V g,. Die Anzahl der Nebenbedingungen vom Typ 9, = 0 
ist also um die Anzahl der Grenzflächen, auf die der Punkt gesetzt worden ist, erhöht worden. 


Ist N die Anzahl der Variablen und M die Anzahl der Funktionen g, und m die Zahl der 
als Nebenbedingungen hinzugenommenen Grenzflächen, so ist zunächst zu prüfen, ob 


N ee EEE 7) 


ist. Ist (17) erfüllt, so kann das Iterationsverfahren zur Erfüllung von g, = 0 durchgeführt 
werden. Ist (17) nicht erfüllt, so ist die Zahl der Variablen zu klein, und die Schrittweite wird 
halbiert. 

Die Ungleichung (17) sei nun erfüllt. Der Punkt r/, der sich jetzt bereits auf einer Grenz- 
fläche befindet, kann dadurch dort festgehalten werden, daß man in dem in den Gleichungen (10) 
bis (12) beschriebenen Iterationsverfahren die entsprechenden Komponenten in den Gradienten 
Null setzt. Damit ist bis zur Beendigung der Iteration gesichert, daß die folgenden Punkte 
und damit auch der Punkt r, auf der Grenzfläche liegen. 


Es ist möglich, daß durch Rechenungenauigkeiten der Punkt r/ nicht exakt bei Auswertung 
von (6) auf die Grenzfläche gelegt wird. In diesem Falle könnte es sich ereignen, daß während 
vieler aufeinanderfolgender Schritte y ständig in ein immer kleineres y* verkleinert wird, ohne 
daß die Grenzfläche wirklich erreicht wird. Um einen derartigen Vorgang zu verhindern, wird 
auch für die Ungleichungen eine Genauigkeit angegeben, die mit b, bezeichnet werde. Ist 


une Bl br Me ee: Fe ern ee 15) 


oder i 
ey Pe A ar A Eee 0 tet:iı) 


so wird angenommen, daß die Koordinate x, auf der betreffenden Grenzfläche liegt. 


Jetzt muß noch der Fall betrachtet werden, daß ein Ausgangspunkt r/, bereits auf einer 
oder auf mehreren Grenzflächen liegt. Es wird dann zunächst wieder entsprechend den Glei- 
chungen (7) und (8) ein Vektor dv, berechnet. Dabei wird jedoch häufig der Fall eintreten, daß 
p, aus dem von den Ungleichungen zugelassenen Bereich dort herausragt, wo der Punkt bereits 
auf der Grenzfläche liegt. 

Es gibt dann zunächst nur die Möglichkeit, die betreffenden Grenzflächen als Neben- 
bedingungen hinzunehmen, das heißt, die betreffenden Komponenten von db, Null zu setzen. 
Damit ist aber im allgemeinen (7) verletzt, da man nicht damit rechnen darf, daß die V 9, ortho- 
gonal zu den Normalvektoren der betreffenden Grenzflächen liegen. Es ist daher notwendig, 
die Auflösung des Gleichungssystems, das sich durch Einsetzen von (9) in (7) ergibt, erneut 
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\ L ;ö u 
a 
vorzunehmen. 'Bezeichnet man mit e, die Einheitsvektoren in Richtung der Normalen 

Nebenbedingung gezählten Grenzflächen, so ist jetzt die Aufgabe zu lösen: 
ee Mr u 
Ra A A en, 
mit 

n=AUVFHZE un V m FZBEe) -» ren 


Wegen (19b) hat d, keine Komponenten in Richtung der e,. Bezeichnet man mit (VP* und # 
(7 9m)* die Gradienten, bei denen die »-ten Komponenten Null gesetzt wurden, so läßt sich 
(20) also schreiben: 


ve A [VDE EVER Re (21). 

In (19) sind damit die v-ten Komponenten ohne Bedeutung, und damit erhält auch (19) die 
Form 

DVG IE. u 00 a a A 


Jetzt ist also das durch (21) und (22) gegebene Gleichungssystem zu lösen, und man erhält 
einen neuen Vektor d,, der möglicherweise bei weiteren Komponenten die Ungleichungen ver- 
letzt, während gleichzeitig die entsprechenden Komponenten von z, auf den Grenzflächen lieg.n. 
Es ist also erneut die Prüfung und erforderlichenfalls die Korrektur von (\Vf)* und (Vg)* in 
((Vf)*)* und ((Vg)*)* vorzunehmen, und das Gleichungssystem (21), (22) ist bezüglich der 
,„ neu aufzulösen. Dieses Verfahren ist solange fortzuführen, bis zu dem Punkt r, ein Punkt r 
entsprechend (6) gefunden werden kann, der keine der Ungleichungen verletzt. Selbstverständ'ich 
ist nach jedem Schritt eine Prüfung gemäß (17) vorzunehmen. Stellt sich heraus, daß keine 
freien Variablen mehr vorhanden sind, so ist die Maximalisierung beendet. 


Gegen das hier beschriebene einfache Verfahren zur Berücksichtigung der Ungleichungen 
läßt sich nun ein Einwand erheben. Es kann der Fall eintreten, daß das Verfahren die Maxi- 
malisierung für beendet erklärt, während noch durchaus eine Möglichkeit bestanden hätte in 
Richtung auf größere Funktionswerte fortzuschreiten. So könnten etwa beim ersten Durchgang 
die Komponenten » zu Null gemacht werden und beim zweiten Durchgang die Komponenten ?, 
und es könnte sich dann herausstellen, daß N— M—»—»=0. Die Maschine registriert 
dann, daß nicht mehr genug freie Variable vorhanden sind und erklärt die Maximalisierung für 
beendet. Möglicherweise ist aber durch die Streichung der Komponenten # die Streichung einer 
der Komponenten » überflüssig geworden, und es bestand daher durchaus noch eine Möglichkeit 
in Richtung auf größere Funktionswerte im erlaubten Bereich fortzuschreiten. 


Dieses Problem soll im nächsten Abschnitt besprochen werden. 


Erwähnt werden soll hier noch kurz die Schwierigkeit, daß durch das Iterationsverfahren 
zur Erfüllung von 9, = 0 die Ungleichungen verletzt werden können. Es kann der Fall ein- 
treten, daß der aus (6) ermittelte Punkt r/ die Ungleichungen (3) noch erfüllt und daß sein Ab- 
stand von der v-ten Grenzfläche auch noch größer als b, ist, während der nach der Iteration 
ermittelte Punkt x, die Ungleichungen (3) verletzt. In diesem Falle empfiehlt es sich, wenn die 
Bedingung (17) es gestattet, die betreffende Grenzfläche als neue Nebenbedingung einzuführen 
oder aber die Schrittweite zu halbieren. Die Iteration zur Erfüllung von 9. = 0 muß dann neu 
durchgeführt werden. 


4. Die strenge Berücksichtigung der Ungleichungen 


Im vorigen Abschnitt wurde ein einfaches Verfahren angegeben, wie die Ungleichungen (3) 
bei der Maximalisierung berücksichtigt werden können. Dieses Verfahren erwies sich insofern 
als unbefriedigend, als bei seiner Verwendung die Suche nach dem Maximum möglicherweise 
vorzeitig aufgegeben wird. Hier soll nun eine strenge Formulierung der Aufgabe besprochen 
und ein Weg zur teilweisen Behebung der Schwierigkeiten gezeigt werden. 


Voraussetzung für das Auftreten der genannten Schwierigkeit war, daß in (6) der Aus- 
gangspunkt r, auf mehreren Grenzflächen lag. Der einfacheren Schreibweise wegen sei ange- 
nommen, daß r, ausschließlich auf Grenzflächen des Typs (3b) liegt. Mit y b soll jetzt der Vektor 


bezeichnet werden, um den x, in den neuen Vektor z, verändert wird, das heißt, an die Stelle 
von (6) tritt jetzt die Gleichung (mit y > 0) 


a a 
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Über den Betrag von vd wird im Augenblick keine Annahme gemacht. Die Gleichungen (7) und 
(8) bekommen jetzt die Form 


BU RAN: 
pn V/= Max RE (05) 


Mit e, seien nun Einheitsvektoren in Richtung der »-ten Koordinatenachse gemeint. Die 
Forderung, daß in den Komponenten », in denen rt, auf den Grenzflächen liegt, der Vektor I, 
die Bedingungen (3b) nicht verletzen darf, läßt sich dann auch formulieren: 


1,020; A re Er PAS) 
(26) muß für alle » erfüllt sein, in denen x, auf einer der Grenzflächen (3b) liegt. 
Es seien nun n, linear unabhängige und orthogonale Einheitsvektoren, die die Gleichungen 
1, V 9,0). =0 Baelessshhe Nele NZ Moe en 
erfüllen. Es gibt N— M derartige Vektoren. Aus (24) folgt, daß sich v darstellen läßt als 
DL N ehe en el 


Es besteht also jetzt die Aufgabe, die a, zu bestimmen, indem man (28) in (25) und (26) einsetzt. 
Sieht man die a, als unabhängige Variable an, so ist also wieder eine Maximalisierungsaufgabe zu 
lösen, die sich von der unter (1) bis (3) gestellten Aufgabe dadurch unterscheidet, daß keine 
Bedingungen der Form g, = 0 vorliegen, daß dafür aber die Ungleichungen eine schwierigere 
Form haben, obgleich sie noch linear sind. 


Die zu lösende Aufgabe hat jetzt folgendes Aussehen: 


I EL er ee euer 
Yza 
ae Er 1; 2, ae ee: 


Auch bei dieser Aufgabe besteht — wenn man sie nach dem in den vorhergehenden Abschnitten 
beschriebenen Verfahren lösen will — die Schwierigkeit, daß möglicherweise die Suche nach 
dem Maximum vorzeitig abgebrochen wird. 

Es gibt allerdings einen Ausweg aus diesen Schwierigkeiten. Dazu muß man bedenken, 
daß es für den Augenblick weniger darauf ankommt, aus (29) ein wirkliches Maximum zu machen 
als festzustellen, ob (29) unter den Nebenbedingungen (30) noch positive Werte annehmen 
kann. Die zweite Frage läßt sich aber auf folgende Aufgabe zurückführen: 


HU RN Er er ee 
Za,(e;n) 20; al Brian Neid) 
na a a Eee tale 


Führt die Aufgabe (31) — die sich nach den Methoden der linearen Planungsrechnung 
bearbeiten läßt — zu keinem positiven Maximum für (31a), so ist mit Sicherheit die allgemeine 
Maximalisierung der Gleichungen (1) bis (3) beendet, da in diesem Fall auch (29) kein positives 
Maximum haben kann. Ist das Maximum positiv, so läßt sich der Wert f noch durch Fort- 
schreiten in der neuen Richtung vd verbessern, wenn auch v im allgemeinen nicht die in (29) 
und (30) geforderte günstigste Richtung sein wird. 


5. Die Benutzung von Ungleichungen allgemeiner Form 


Wesentlich allgemeiner wird die Aufgabenstellung, wenn man an Stelle der Ungleichungen 
(3) Bedingungen von der Form 
een ner) 


zuläßt. Zunächst ist festzustellen, daß die Anzahl der Bedingungen der Form (32) theoretisch 
keinerlei Beschränkungen unterliegt. Man denke als Beispiel daran, daß das Maximum einer 
Funktion im Innern eines auf einer Ebene gelegenen Vielecks zu suchen sei. Dann besteht pro 
Seite des Vielecks eine Ungleichung der Form (32). 


u 
ME “ ei n 
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Nimmt man eine Erhöhung der Zahl der Variablen in Kauf, so lassen sich ( 
gemeinen Ungleichungen durch Gleichungen der Form (2) und Ungleichungen der Fo N 
ausdrücken, wenn man pro Ungleichung der Form (32) eine neue Variable xy,, einführt: 


i Rh, ar; ed, 
ml: 


Will man keine zusätzlichen Variablen einführen, so kann man auch mit den Ungle’chungen (32) 
direkt rechnen. nr 

Ausgangspunkt sei wieder ein Punkt im Innern des erlaubten Bereichs. Nach dem bereits 
beschriebenen Verfahren kann man dann solange in Richtung des Gradienten fortschreiten, 
bis man auf eine Grenzfläche stößt. Ist h, keine lineare Funktion, so muß beim letzten Teil 
der Annäherung an die Grenzfläche ein Iterationsverfahren verwandt werden. 

Beim nächsten Schritt der Maximalisierung befindet sich nun der Ausgangspunkt bereits 
auf der Grenzfläche. Weist der Vektor v, dann bezüglich dieser Grenzfläche in den verbotenen 
Bereich, so ist zunächst wieder eine Prüfung analog (17) auf die Anzahl der noch zur Verfügung 
stehenden Variablen durchzuführen, und anschließend nimmt man die betreffende Grenzfläche — 
wenn die Prüfung ergab, daß noch freie Variable vorhanden sind — als Nebenbedingung der 


Form (2) hinzu. Das geschieht dadurch, daß man in dem Gleichungssystem (19a), (19b) und u 
(20) den Vektor e, durch den Gradienten von h, ersetzt. Eine Vereinfachung, wie sie durch den £ 
Übergang auf (21) und (22) eingeführt werden konnte, ist jetzt nicht mehr möglich. ar 


Liegt bei späteren Maximalisierungsschritten der Ausgangspunkt auf mehreren Grenz- 
_ flächen, so ist analog zu verfahren. 

Ergibt bei einem der späteren Fälle Gleichung (17), daß die Zahl der Variablen nicht mehr 
ausreicht, so ist der in den Gleichungen (23) bis (31) beschriebene Weg einzuschlagen. Bezeichnet 
man wieder mit » die Indices der Grenzflächen, auf denen der Ausgangspunkt bereits liegt, so 
tritt an die Stelle der Vektoren e,; in den Gleichungen (23)—(31) der Gradient der Funktionen A;. 
Man findet auch jetzt wieder eine günstige Richtung zur weiteren Maximalisierung nach einem 
Verfahren der linearen Planungsrechnung, und man erhält eine zutreffende Antwort auf die 
Frage, ob die Zahl der Variablen zu einer weiteren Maximalisierung noch ausreicht oder nicht. 


6. Bisherige Erfahrungen 


Es wurde für die IBM 650 ein Programm aufgestellt, das den Teil des Maximalisierungs- 
verfahrens enthält, der in den Abschnitten 1 bis 3 beschrieben ist. Das aufgestellte Programm 
gilt fürN= Wund M=8. 


Bei den Fällen, in denen das Programm bisher angewandt wurde, hat es sich gut bewährt. 


.7. Zusammenfassung 


Es wird ein auf die Praxis zugeschnittenes Verfahren beschrieben, wie man sich dem 
Maximum einer Funktion annähern kann, für deren Variable Nebenbedingungen in Form von 
Gleichungen und Ungleichungen bestehen. Das Wesen des Verfahrens besteht darin, daß man 
einen vorgegebenen Punkt in systematischer Weise so verschiebt, daß man ohne Verletzung 
der Nebenbedingungen zu immer höheren Funktionswerten gelangt. Irgendeine Aussage über 
den Grad der Annäherung an das Maximum kann zu keiner Zeit im Rahmen dieses Verfahrens 
gemacht werden, solange man nicht über die benutzten Funktionen neben der Differenzierbarkeit 
weitere Voraussetzungen einführt. 


Ernste Schwierigkeiten treten bei dem Verfahren erst durch die Berücksichtigung von 
Ungleichungen als Nebenbedingungen auf, und zwar frühestens, wenn man den Punkt auf 
seinem Weg zum Maximum auf mehr als eine Grenzfläche gleichzeitig legt. Der für diesen Fall 
gezeigte Lösungsweg ist noch nicht befriedigend. Eine einfache Form des beschriebenen Ver- 
fahrens — einfach vor allem in bezug auf die Art der zugelassenen Ungleichungen — wurde 
für die IBM 650 programmiert und hat sich in einer Anzahl von Fällen praktisch gut bewährt. 

Ich danke Herrn Dipl.-Math. G. Viebeck von der mathematischen Abteilung der IBM in 
Sindelfingen und meinem Kollegen, Herrn Dipl.-Math. W. Ziegler, für die gute Zusammen- 
arbeit auf dem Weg von der ersten, ungenauen Formulierung bis zum fertigen Maschinenpro- 
gramm. Mr. D. Walley von Babcock and Wilcox Ltd., London, danke ich für anregende Rat- 
schläge. Der IBM-Deutschland danke ich für die freundlicherweise bereitgestellte Maschinenzeit. 


Manuskripteingang: 19. 10. 1958 


Anschrift: Dr. P. St. Pütter, i. Hs. Deutsche Babeock & Wilcox-Dampfkesselwerke AG, 
Oberhausen (Rheinl.). Postfach 34/35 
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Unsteady motion of a viscous liquid contained between two 
Bir: infinite coaxial cylinders r=a and r=b 


By C.D. Ghildyal 


Für eine zähe homogene inkompressible Flüssigleit, die sich zwischen zwei koaxialen Zylindern befindet, 
wird in der vorliegenden Arbeit eine nicht-stationäre Bewegung gefunden. Das Resultat stellt eine Erwei- 
terung des von Sneddon (1951) angegebenen Ergebnisses dar. Anfänglich rotiert der äußere Zylinder mit 
konstanter Winkelgeschwindigkeit 2. Wenn der stationäre Zustand erreicht ist, wird der innere Zylinder 
spontan in Drehung versetzt, so daß er mit konstanter Winkelgeschwindigkeit Qı zu rotieren beginnt. 


In this paper we have obtained an unsteady motion of a viscous homogeneous incompressible liquid, 
contained between two coaxial cylinders. The result is an extension of that given by, Sneddon (1951). Ini- 
tially, the outer cylinder is rotating with constant angular velocity Q. When the steady state is attained, 
the inner cylinder is given an impulsive twist such that it begins to rotate with constant angular velocity Aı. 


Dans l’article present on trouve un mouvement non-stationnaire pour un liquide visqueux, homogene et 
incompressible, se trouvant entre deux cylindres coaxiaux. Le resultat represente une extension du resultat 
indique par Sneddon (1951). Au commencement le cylindre exterieure se trouve en rotation avec une vi- 
lesse angulaire constante 2. Au moment ou l’etat stationnaire est atteint le cylindre interieure recoit une impul- 
sion rotative, de sorte qu’il commence a tourner d’une vitesse angulaire constante Q4. . 


B HacToamei pa6oTe HolyyaeTcA HeyCTaHoBHUBlIeecH HBH;KeHNE BA3BKOÜ, OMHOPONHOH 
MH HEeCHKHMAaeMoÜ 3KHIKOCTH, COMep;KalleficH MeKNy ABYMA KOAKCHAJIbHbLIMU HMAIHHNPAaMuH. 
Pe3yJIbBTaT ABAIAETCA 000Ö0INEHHEM pe3yIbTaTa, NOIy4YeHHOTO (B 1951 Tr.) CHejtoHom. B Hayalıe 
BHEINHHÜ IMAHHNP Bpamaertca C IIOCTOAHHON YTIIOBOAÜ CKOPOCTBIO 2. ITlocne NOCTWKEHUA 
YCTAHOBHBIHETOCH COCTOAHUA, BHYTPeHHeMY HHJIMHAPY NaeTcA BpamaTejıbHbIH TOIUOK, TAK 
yYTO OH TaK}Ke HAYHHACT BPAIIATbCA C HOCTOAHHOH YTIIOBOÜ CKOPOCTBP 21. 


Introduetion 


The steady motion of a viscous homogeneous incompressible liquid between two coaxial 
cylinders has been discussed by various authors, but only a few unsteady solutions are known at 
present. To solve the equations of motion in unsteady cases, it was felt necessary to make use of 
the various known transforms. 

Sneddon (1951) in his book ‚‚Fourier Transforms‘ discussed the application of some well 
known transforms to the problems of applied mathematics. He also obtained an unsteady solu- 
tion when the outer cylinder is rotating with constant angular velocity 2, and the inner cylinder 
is kept fixed. Mallick (1957) studied the unsteady case when an infinite circular cylinder is given 
an impulsive twist, at time £ = 0, such that it starts rotating in an infinite, viscous and steady 
rotating liquid, with constant angular velocity 2,. The axis of the rotating liquid and the cylinder 
is the same. ’ 


1. Formulation of the Problem 


Let the viscous homogeneous incompressible liquid be contained between two infinite 
coaxial cylinderssr=aandr=b(b > a). Let the outer cylinder be rotating with constant an- 
gular velocity 2 and the inner cylinder be fixed. In cylindrical polar coordinates (r, ©, Z) the 
equations of motion reduce to 


LE: 
TE 0.68 Or, 
and 
ov dv 1 v 
Er RU er nn P (10%: 
ot dr Ta ) 0) 
In the steady case 
Dv= Br 
Evaluating A and B, with boundary conditions 
erirrearand Wo=ßb, r=b, \ 
we get Te i 
Barcee Fa 


r &—a)' 


Let us suppose now that an impulsive twist is given to the inner cylinder so that it rotates 
with constant angular velocity 2:- 
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The equation of motion then is 


ov dal Ed 
—_ = —— ——.,... 2.2.2... ren 1.2 
ot dr r or | (1.2) 
with boundary conditions 
() Bann © 
(ii) v=bQ, r=2 (1.3 
a are N late ee N}: 
Tr (#®—a)’ En 
2. Solution R 
Multiplying equations (1.2) by!) 
rl) Yıbd 5) — Yılr 5) Jılb &)} 
where £, is a positive root of the equation 
Klac)r &EETBE Vıla a) Or (2.0), 
and integrating from a to b, we get 
dv Wine Jı(b &;) = 
a aan 0 \ nk Den er (2.1) 
where 


[) 
u= e% riIr Ei) Yılb &) — Ib 5) Yılr5)} var. 


Integrating (2.1), we get 


AR en. 
D = & Er? Et 2 +ce 
also 
= 2b 
D Iı =0o  —— & . 
Therefore, 
ES BAZD) 
RT 
Hence 


Ev =—bQ + el EA 
1 i 


By using inversion formula, we get 


r BB — a r DB _—_ (a 


le ee ee 
IE) ag) FH YOH- er Ab. 


D 


2a, 


Using (2.0), this can be written as 
_2u#r—a D,0(b?—r) 


v 22? 
r 2b — a2 r (b? — a?) 


2 BE) ei’ : 
+2a0, 2 A 2) TE a) Vila = I+AHT.. (22), 


The value of v satisfies the boundary conditions (1.3) and also the equations (1.2). 


‘) Jı(@) and Yı(z) are Bessel Functions of order one of the first and second kinds, respectively. 


# 
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I is the initial velocity and A is the part added to it when the steady state is finally attained. 
T is the transient part which ultimately disappears. The behaviour of T for any fixed ratio — ; 


say k, can be easily discussed. Let a be the unit of length, then b=k. The value of T now 
becomes 


1=29,3 ERDE" re) EL 
Ne) J2(k €) i ı\si i) Jı\Si 


It is now easy to see that 


(1) if the value of k is quite near to unity, the steady state is attained very soon after the disturb- 
ance. 


(ii) if k is very large the unsteady state will continue for a longer period. 


For values of k other than those considered above, let 7 be the smallest root of (2.0), then 
we have 


— re Jik E) {Jr &) Yıld) — Yılr &) JE} MER anı 
190% a 11 Si) ü ı\si i) 41\Si + e 
re 2, RE) — RE) - a 
Therefore at any fixed time i,, T decreases as r increases from Are to k(k >.1). 
The velocity at any instant Z and at distance r from the axis of rotation is given by the 
expression (2.2). 
If, in particular, we take »t = 0.1 and k = 2, we get 
ET 
e 1 —- = —_ 
272 (3-19) — I? (6 - 38) 
OS DL Te le 
7 (6:30 72.0) 
+ +. + smaller terms. 
Ignoring the smaller terms and simplifying the above expression, we get 
T = K, {3688 J, @ - 19 r) — 2654 Y, 8 -19r) 
+01 x Y,(6-317) — 0102 x J,(6-31r)} +: --, 
where K, = 0.704. 2.. 


We give below the values of T for different values of r(1<r<<2). Itis easy to see from 
the table that T decreases from zero atr=1,toa minimumatr=r,(r,< 3/2). For, <r<2, 
T increases steadily to zero. 


122.0 {J,@-19r) Y,@-19)— Y, (8-19 r) J, 8 - 19)} 


{Jı (6-31) Y, (6-31) — Y, (6-31) J,(6-31)} + - + 


Table showing relation between r and 7 


7 | 


r T r 7% r | Ih 


| 
| I 
1 .0.00000 47/322 | —0.16831 K, | 31/16 | 0.027339 X, 
33/322 | —0.02083 K, 3/2 —0.16695 K, | 63/32 —0.01345 K, 
17/16  —-0.04068K, 49/32 —0.16351 K, | 127/64 —0.00661 K, 
98 | —0.07611 KL 25/16 | —0.16127 K, | 255/128 | —-0.00388 K, 
54 | —0.13294 K, | 13/8 —0.14544 K, ı 511/256 —0.00108 K, 
1150 00.163534. K, 7/4 —0.13439 K, 0.00000 
23/16 | —0.17074K, 15/8 —0.05444 K, | 


3. The friction couple per unit length of the cylinder r = a, is given by 


where 


The expression (2.2) can also be written 
„al 8 a0 (0, — 2) 
B—a r (b? — a2) 


2a an {Jı(r &) Yıla &) — Jıla &) Yılr &)} - 
US 9 (a E) J2(b £,) Si i 1 i 1 i 
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This gives 


8. fv} 
se er 7 


& 20289, —Q 200, zone“ 
u a a ee 1 


6 Ra Yıla) Aa 5) Yılr 50) 


209 m Met. Br | an 
HG Fa) RE Ya) — Ye) ns]. 


Hence 


20,9) ZI 
Prol=a = 24 I- "(be eu +2 2 Ja &) —Ji(b 5) 


Thus we get?) 


Met none 
G=4anu |- per +2, 2 JS: (a&)— J?(b | 


In the case when a = 1,b = 2andvit = .1, we get 


(4 tl, 9 + 37x 2 


approximately. 


I am grateful to Dr. Ram Ballabh for supervising this work and to Scientific Research 
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Ein Beitrag zur Dynamik der Kreismembran 
Von Max Hieke 


Mit Verwendung des Dirichletschen diskontinwierlichen Faktors wird eine inhomogene Differential- 
gleichung für die Auslenkung der Kreismembran integriert. Dabei wird vorausgesetzt, daß die Störungs- 
funktion auf einem Segment der Kreisfläche konstant ist, aber auf der restlichen Kreisfläche verschwindet. 

Dirichlet’s discontinuous factor is used to integrate the inhomogeneous differential equation for the 
displacement of a circular membrane. T'he disturbance is supposed to be constant over a segment of the mem- 
brane and to vanish elsewhere. 

En appliquant le facteur discontinu de Dirichlet une equation differentielle non-homogene pour le de- 
placement de la membrane circulaire est integree. On y suppose que sur un segment de l’aire du cercle la 
fonction des perturbations soit constante, tandis qu’elle disparaisse sur le reste de l’aire du cercle. 


IIpn nomoımm IPepbIBHCTOrO MHOSKMTENA JInpuxie HHTETPuUPpyeTcH HeomHoponHoe nuhehe- 
PeHnmasIbHoe yPaBHeHNHe UA OTKIOHCHNA KPYTOBOH MeMÖpanpı. Ilpı 9ToMm npennonaraetcn, 


UTO $bykımsı BOBMYINeHMI HA OJAHOM CETMEHTE KPYTAa IIOCTOAHHA, aHa OCTAalIbHOH yacTu Kpyra 
paBHa HyJIIO. 


Für die Untersuchung bestimmter Schwingungsformen einer am Rande fest eingespannten, 
kreisförmig berandeten Membran mit dem Radius r = a sei der Mittelpunkt 0 des Randkreises 
der Mittelebene der unbelasteten Membran zugleich Ursprung eines ebenen Polarkoordinaten- 
systems {0; r, p}. Die vertikale Verschiebung w(r, 9, t) eines beliebigen Punktes der Membran- 
mittelebene wird unter transversaler Belastung durch die inhomogene partielle Differential- 
gleichung 
1 &ofr, 9, t) 


Aofr, 9, 1) — _ 
on pd)— Er 


= Bor, DH 2A De 
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mit der Randbedingung 
HT, E en e  R, 20, Wr aiee 


bestimmt. Darin ist B= — > gesetzt worden, und es bedeuten weiter c = /T/« die Schall- 
42 ’ 


ausbreitungsgeschwindigkeit innerhalb des Membranmediums und 4 = «ad die Membranmasse 
je Flächeneinheit. T ist die mechanische Spannung, d die Dicke und & die Dichte der Membran. 
Die Funktion p(r, 9, f) ist die auf die Flächeneinheit bezogene transversale Belastung der Mem- 
bran. 


Es sei eine zur Zeit £ = 0 plötzlich einsetzende Belastung der Membranfläche angenommen, 
dargestellt durch die Funktion 


Pı(r, 0), fürate>0 


a er Se a ET, 


iD) 
Als Lastverteilungsfunktion p,(r, @) sei die gleichförmige Belastung eines Kreisabschnittes 
des Randkreises der Membranoberfläche gewählt. 


NP für rcsp—1>0 
Al 9) = Nie für rcosp—I<0O. 


Die spezifische Belastung ist danach nur innerhalb 
eines Kreisabschnittes vom Öffnungswinkel 29, der durch 
r = a begrenzten Kreisfläche von Null verschieden, wäh- 
rend sie im übrigen Teil der Fläche verschwindet (s. Bild 1). 
Für den Winkel ©, gilt die Bestimmungsgleichung 
acsay —I1=0. 

Mit Hilfe des Dirichletschen diskontinuierlichen Fak- 
tors läßt sich die vorgegebene Belastungsfunktion p(r, 9, t) 
wie folgt darstellen 


B Po [* eitt de eiz(r cosp—l) 


p(r, 9 t) = An? 


dz (1b), 


Bild 1 


ce 


wobei die Integrationswege c bzw. c’ je längs einer Parallelen zur reellen Achse der komplexen 
z- bzw. ö-Ebene unter Vermeidung der singulären Stellen zu erstrecken sind. 


Der inhomogenen partiellen Differentialgleichung (1) in Verbindung mit der Randbedin- 
gung (la) genügt formal der Ansatz 


B zeit? e-izl = +0 dz 
oT, 9, f) ae © de | RE ee EN (KT) el ee (2), 


v=—oo 
et 


worin abkürzend gesetzt ist 
2# 


J eiazcosd end d®. 


S 2 Bari 
2 Ilka) , 


k= a. = 
7 
G 


Dabei ist der erste Summand dieses Ansatzes (2) ein Integral der inhomogenen Differential- 
gleichung (1), während die unendliche Reihe formal der zu (1) gehörigen homogenen Differential- 
gleichung genügt. Setzen wir voraus, daß sich die Reihenkoeffizienten so bestimmen lassen, 
daß diese unendliche Reihe gleichmäßig konvergiert, so darf unter Vertauschung der Integrations- 
reihenfolge gliedweise integriert werden. Unter dieser Annahme finden wir aus dem inneren 
Linienintegral von (2) mit der Definition des diskontinuierlichen Faktors 


[1 für resp —I>0 


== 2.08 
& \o für rcsp—1<0 @) 
den Wert 
2ril= |. © 2nilat1-cosk(rcosp—D}E 
er +® 
1 X Jlkn air 24 —ivd 
1 TE. TR‘ p | {—1+cosk(acos#®—D)}e dd|. 


—ps 
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Mit dem vorstehenden Ausdruck erhalten wir für die Membranauslenkung zunächst das 


rail 
Linienin egr Eu f (k = &jo) et ” 
a an 


ET 2rni 


ä i -Ebeı ter Vermeidung der sin- 
das längs einer Parallelen zur reellen Achse der komplexen Z-Ebene unter 
sulkten Stellen des Integranden zu führen war. Da die Funktion I) sich asymptotisch dem 
Wert Null nähert, wenn £ in seiner positiv-imaginären Halbebene über alle Grenzen wächst, 
sind die Voraussetzungen für die Anwendbarkeit des Residuumsatzes auf das vorstehende 
Linienintegral erfüllt. \ 

Um die Membranauslenkung w(r, 9, l) zu gewinnen, berechnen wir zunächst das Residuum 
an der Stelle &= 0. Wir finden dafür 


+ 


’ 1 nn Ir > iv 
app a ut x; e j (a cos d — I)? e 9) dd (4). 


2ri an, a 


=—o 


Po 


Die erhaltene Reihe läßt sich reell schreiben und stellt die im Innern des Kreises r <a 
gültige Fourierentwicklung des längs des Randkreises r = a erstreckten Poissonschen Integrales 
dar, das über den Randwert (acos#—)?& zu nehmen ist. Denn man hat 


vw=+® 
1 


+9 
Pa A 0 .a—r (acosd— dd 
2n 2 (2) | er le er I u 


v=—o 


— 9 Ps 


Die Auswertung des Poissonschen Integrales ergibt die folgende Funktion 


+ 


a? — r2 (a cos d — I)? dd 2 
In j a2 1 r—2arcos(—9) 


Ps 


DE hr RT 


mit der folgenden Festsetzung des diskontinuierlichen Faktors 


0 für —_n <p <— 
CHstl für m RL real a ae 
0 fü +9, <po<+n 
In (5) ist nun 


2 ___ r2\2 2 2 
D— Kanes cos29 + aco?p +1 — lcosp a: ir) 
4r2 r 
zu setzen, und die Funktionen ®; = Ö(r, 9, d;) miti=1,2, und , = 90 = — 9, ‚lauten 
a — r? at — rt 
) = — |a? cos 2 9 1— j —d) + — —9d; 
Ö(r, o, d;) a7 |@ cos2 NE sin (9 — 9;) 4 mn (® — d;) 
(a? — r2)? d—r 0 a—r en 
u Yan arcig in: cot 5 + 2 (a? cos?p +12) arc lan „ct 3 Pr 
z ge a? — -ri A A iz e r | at 3 ri | ’ 
+2alsinp m In RE —a?sin 2 TE cos (P— d)+ dan In R; 
a — r? s a@—+r kr re 
-2aleosp|,_. —9d) + m arc tg ne p) ) 
wobei 


R=a+nr—2arcos(p— 9) 
zu setzen ist. 


Damit erhalten wir für den statischen Anteil der Durchbiegung der Membran, den wir der 
Kürze halber mit wu(r, p) bezeichnen, den Ausdruck 


B 
a) = OP MC-BC—-DHD........0, 


Darin hat & die nach der Festsetzung (3) und &’ 
Durchbiegung (Tr, @) ist eine im Innern r < a 
Funktion, die in Erfüllung der geforderten Ran 


die nach (6) gegebene Bedeutung. Die statische 
und auf dem Rande des Kreises r — a stetige 
dbedingung für r = a verschwindet. Identifi- 


. 
N \ CE. vn 
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ziert man in diesem Ausdruck I! mit (— a), so erhält man die Durchbiegung der über ihre gesamte 
Fläche gleichförmig belasteten Membran. Jetzt sind noch die den Residuen an den einfachen 
von Null verschiedenen Wurzeln der Besselschen Funktionen entsprechenden Anteile des 
komplexen Integrals zu ermitteln. Für diese Nullstellen gelten die Beziehungen 


t Jn(kum a) = JIn(Oam) = 0, 
wobei 
Snm a 


= EOnm» Khmt=Om= al. 


die Wurzeln der transzendenten Gleichungen sind. Insbesondere setzen wir für die Besselsche 
Funktion nullter Ordnung J,(o) 


Oom = Om» Om = km a = —. 


Um eine handliche Bezeichnung zu haben, sei der erwartete Ausdruck mit W bezeichnet. 
Man erhält dafür die Doppelreihe 


BEN Ne De rote 
RT 0 a em 
RA) 


Da der von den Nullstellen © = 0 der Besselschen Funktionen herrührende Anteil bereits 
ermittelt ist, ist bei der Summation über u der Wert u = 0 auszuschließen. Diese Vorschrift 
ist durch den Stern am Summenzeichen ausgedrückt. Die Konstanten A,, haben den Wert 


A, a ee 


v=—o u=—o 237 


Ip 


+% 
Ayn = {—1+cosk,,(acs# —D)}e"?’dd = = sin» gt —1-+ cos (k, „DJ. a} 


PP 


= 2 cos (Ku D) 2 eo Dr JIorlkyu a) 


sin@r+)M , sin@r—») & (11). 
2T7+» 2T—v 

sin @r Fl 69) 9% , sin (2r + I 2. 
2r+1+v 2r+1—» 


+ 2 sin (k,„D) >, (— 1)" Jarzılk,. d) 


Die Doppelreihe (10) schreibt sich reell so: 


BP a? S Jolku r) Ay T JIckyu r) Ay 
W= —- 7 6o5L,E+2 ——— —-cosvopcosl,„ti. . (10). 
7 2, | 04 J6(0,) 3 2 Or u Jr(0, 1) : 


Die Reihendarstellung der Konstanten A,„, entsprechend » = (0, erhält man aus (11) 
durch Ausführung des Grenzüberganges für v—0. Die Funktion W(r, 9, f) stellt die freien 
Schwingungen der auf einem Kreisabschnitt statisch angeregten Membran mit fest eingespann- 
tem Rand dar. Die gesamte Durchbiegung der Membran wird durch die Funktion 


&1,.0D =@lt.0) + Wird... +. 2 Meer 7 2) 


beschrieben. Man erhält die Reihenentwicklung der Funktion — o,(T, 9), mit w,(r, 9) als dem 
statischen Anteil der Ausbiegung, wenn man in der Funktion W(r, o, t) die Zeit {= 0 setzt. 
Es ist 

{W(r, 9, Dh=o = — w(T, P) - 


I. Unterfall 


Die Ausbiegung einer gleichmäßig über einen parallelen Streifen belasteten Membran 
erhält man aus den obigen Formeln, indem man den Parallelstreifen als Differenz zweier Kreis- 
abschnitte auffaßt und mit den dazu gehörigen Ausdrücken (9) und (10) für die jeweilige Aus- 
biegung entsprechend verfährt. Ersetzt man den gemeinsamen konstanten Faktor p, durch den 


Quotienten 
Po = PılAl 


und geht in der Differenz beider Funktionen w(r, @, f) zur Grenze über, indem man mit Al ZUu- 
gleich die Streifenbreite nach Null gehen läßt, so erhält man beim Grenzübergang aus dem Dif- 
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ferenzenquotienten beider Ausbiegungen w(r, 9, t) die zur Linienlast gehörende Ausbiegung 
&(r, 9,1). Die Ausbiegung bei Linienlast ergibt sich danach durch Differentiation der Funktion 
o(r, 9, t) (12) nach dem Abstand !. Bildet man zunächst die Ableitung von der Funktion (9), 
so erhält man für den statischen Anteil der Durchbiegung bei Linienlast die F unktion 


Be BB NT 
IR) 5, t2resp—_DE-XE — N, ER, T IIMEIERE BE 
Darin sind 
zung a — ee 
X; In 4 larctg ee‘ 9 +2asing re In Ri 
ner a Eh 
—2ac0sg| Sr —d) + 5 arctg ee | .. . (14a) | 
als Ableitung des Ausdrucks (8), und 
2 z 
X=2i- umge in unlzel 3, rl, VE 


als Ableitung des Ausdrucks von (7) zu setzen. Die hierbei eingeführten oder wieder verwendeten 
Funktionen haben die folgenden Bedeutungen 


X; —_ X,;(r, ® d;) 
mit 

R=a@+nr-—2arcos(p —d,) 
und d, = 9 da = — Po- hr 

Auch der dynamische Anteil der Lösung bei Linienbelastung ist aus der Lösung für den 

gleichförmig belasteten Kreisabschnitt ableitbar. Man erhält ihn aus dem dynamischen Anteil 
W der Lösung (12) bei gleichförmig belastetem Kreisabschnitt durch Differentiation nach dem 
Abstand I. Da! lediglich in den Koeffizienten A, und A,, erscheint, genügt hierfür die Berech- 
nung der Differentialquotienten für die Koeffizienten (11). Man erhält durch Bildung der Ab- 
leitung von (11) für die A,, den Reihenausdruck 


— oA, sin 9 : 
A, ee ku Sin (kun D Jo(kyu 0) 
— 2 ke Gin (ep) uk 1) Jrafk,, a Er Se EEE 
z=1 2trt+rv 27 —v 


sin@r+1+»)p , sin (27 +.1—9)9 


2r+1+4v 2 +1— (0; 


+ 2k,, 605 (k,„) 3 (— VD" Jarrılk,„a) 
z=0 


Diese Darstellung umfaßt zugleich diejenige der Koeffizienten AsR = u Man erhält letztere 
aus der allgemeineren Darstellung (15) der A,, durch die Ausführung des Grenzüberganges für 
v—0. Anstatt des Faktors p, in W(r, 9, !) ist im daraus abgeleiteten Anteil W(r, 9, f) die Linien- 
last P, zu setzen. 

Möchte man den Sonderfall einer längs eines ihrer Durchmesser gleichmäßig belasteten 


Kreismembran erhalten, so hat man in den vorhergehenden Formeln 9, = ” und 1=0 zu 
2 


setzen. Dadurch erzielt man aus (15) für die Reihenkoeffizienten den Ausdruck 


Mei rar ni 
u vu cos v 32 Jerzılk,u a) (27 a” 1% —y3 Eier A (16). 


Bei dieser Anregungsart verschwinden alle Oberschwingungen von ungeradzahligen Index »: 
denn für alev=217+1,7’=0,1,2,... sind die Koeffizienten ärg gleich Null. 


Die längs eines Durchmessers plötzlich durch eine Lini ) i 
| ) ss enlast erregte Kreismembran mi 
eingespanntem Rand führt nur solche Oberschwingungen aus, bei ES; der angeregte RB. 
messer nicht Knotenlinie ist. Aus (13) erhält man die dazugehörige statische Auslenkung 


= BP, 
oy(T,p) = z {2rcsg&— X*Ü"—- XI + X)... 2.7, SH173: 
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mit den folgenden Definitionen 


X 2 2ame[, Be 94T aretg(? m ) 
Bern en ae 
2ar 
X*= — ac0osp —— Dan ee ein 
Wr Di area 2 


Der dynamische Anteil der Durchbiegung lautet 


BP,a Sı [Julkur) A, S J(kuur) Ar 
W'(r, 9,1) = = je - cost,t+2 > -— cosvo cost,,t 19). 
N ren 


II. Unterfall 


In allen bisherigen Untersuchungen war die Polarachse des benutzten Koordinatensystems 
zugleich die Symmetrieachse der Lastverteilungsfläche der Kreismembran. Wegen der Sym- 
metrie der belasteten Gebiete konnte das ohne Beschränkung der Allgemeinheit geschehen. 


Lassen wir die Symmetrieachse nicht mit der 
Polarachse zusammenfallen, so haben wir in den 
Formeln (7), (8), (9) u. (10) anstatt des Winkels o die 
Winkeldifferenz  — & zu verwenden. Der Winkel a 
wird von der Symmetrieachse und der Polarachse 
eingeschlossen. Andere Änderungen der ursprüng- 
lichen Ausdrücke ergeben sich nicht. Die bisherigen 
Ergebnisse lassen sich durch diese einfache Vorschrift 
dahin verallgemeinern. 


Subtrahieren wir zwei unter verschieden ge- 
richteten Symmetrieachsen @ = &;;i = 1, 2 liegende 
inhaltsgleiche Kreissegmente voneinander, so erhalten 
wir als Restfläche zwei Kreisdreiecke, deren Flächen-- 
inhalte gleichen absoluten Betrag haben. Die beiden 
Dreiecke liegen spiegelbildlich zur Halbgeraden 9 = 


= (0 — %,), auf der zugleich ihre gemeinsame Spitze 


Bild 2 


liegt. Ihre Umlaufsrichtungen sind entgegengesetzt 
orientiert (s. Bild 2). 


Lassen wir auf beiden Kreisabschnitten die gleiche konstante spezifische Flächenlast p, 
wirken, so fölgt aus dem Superpositionsprinzip die Lastverteilungsfunktion einer Membran, die 
auf den zwei Dreiecken entgegengesetzt gleich beansprucht ist. Dieses Kräftepaar erzeugt eine 
Ausbiegung der Membran, die gleich der Differenz der Ausbiegungen ist, die zu den einzelnen 
Lastverteilungen über die Segmente gehören. 


Ersetzen wir in der Differenzfunktion die Flächenbelastung durch den Quotienten 


Po asin 9 Aa 

mit P, als Linienlast, und bilden den Grenzwert dieses Differenzenquotienten für da =, —,—0 
so erhalten wir die Durchbiegung einer Kreismembran, die längs einer Sehne so beansprucht 
wird, daß deren eine Hälfte gleichmäßig positiv und deren andere gleichmäßig negativ belastet 
ist. Die Lastverteilungsfunktion bildet ein linienhaftes Moment. 


Den zu dieser Belastungsfunktion gehörenden statischen Anteil der Durchbiegung für die 
am Rande eingespannte Membran erhält man aus der Funktion (9), wenn man darin die Winkel- 
variable og durch  — a ersetzt und den so erhaltenen Ausdruck nach « ableitet. Es ergibt sich 


dann 


—— [2 rsin (@— 0) {res —o)— 1} — 8,® — P9,,+ PD.) (20) 


Boa (T, 9,6) = 2 asin en 
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mit den aus den Relationen (7) und Mi folgenden Definitionen 


Bu=—.,- 2asin2p—) [5 Ar N: ey: a # yı) 
+ ua tg Fe cot en 
—20c0s2 (pP —o) ar ent Tan" Ri) 2 
+ 2alsin Pa) SYr u sr "FE aretg (1 ent? in na 
—2alcos(—o) Fu le RI—2 0° sin 2 (p—a)arctg(“ Sr Zen: 
o, an) nn alsin (pP —o) + @sin2 (p— a) 


Die diskontinuierlichen Faktoren € und 6’ sind nach den Definitionen (3) und (6) zu verwenden. 


Die dazugehörige dynamische Auslenkung folgt aus der Beziehung (10), wenn darin die 
Winkelvariable durch — a ersetzt und der so entstandene Ausdruck nach «& differenziert 
wird. Man erhält dann 


%& 


BP,a | = J,(k,.r)» | 
ee 2 de A,,sinvy(g—o)cost,,t - - - - . (22). 
an % u J»(0, „) i v i 


III. Unterfall 


Überlagert man diesem Belastungsfall den Fall der längs der gleichen Kreissehne kon- 
stanten Linienlast vom Betrage P,, so erhält man den Belastungsfall einer nur längs der halben 
Sehne mit der Linienlast P, = 2 P, angeregten Membran. Insbesondere erhält man für g, = x/2 
die Verhältnisse einer längs eines ihrer Radien angeregten Kreismembran. 

5 auch ! = 0 gilt. 

Die dynamische Auslenkung mit der Polarachse als Symmetrieachse ist für diesen Fall 

nach (19) u. (22) 


BP a S[J,(k,r) A, 
W* + Wi, = > kur) cosö,t+ Re: (ku) Ar "cosvgeosz,,il 


Hier werden auch die Reihenkoeffizienten einfacher, da mit @, = 


| Je) — Bulle. 
BPa<x = v J(k,ur) A 
+ — u 21 ram v m COS il ER a DE Pe 
“7 2 =; Or EEE z vn 22 


mit 


£ PART. = v i 
4, sinv {14 Joh} + 42, DH Il dar (2) 


und der vorstehenden Bedeutung (16) der Reihenkoeffizienten A, für den Fall II. Die zu der 
jetzigen Belastung gehörige statische Durchbiegung lautet nach der Zusammenfassung der ent- 
sprechenden Ausdrücke (17) und (20) 
5 BT, BA 
=, + wo, a0 = 3 {Qarcosp+ rRsn29C— HU’ —H, + Hr >. 99 
mit den folgenden Definitionen 


Es a+r 


1 a—r 
9 El p- dv) + Bu er T eot ?— =) +2@sinp Sn! 


s a? — F 2)2 Be Fa, 
+2@sin2p| 9a B sin) +T - + en arctg(" „car ?Z%) 


a-+r 2 


n a—r } 
2 a®2sin2 t ( ) 
"Fr ınZ2garctg EL „eos 9 » 


dar 


2_r 
+2a@cos2g]% "005 (p— en 


9a =, sei R} 


4a I. 


a? T® a — 2\2 
H=—arcosp“ + @sin 29 57) + asin29 
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und deren Abkürzungen 


R=@+nr—2arcos(— 9); me »=—5- 


Mit der Bestimmung der Lösung der Schwingungsaufgabe für eine Kreismembran, die auf 
einem Segment gleichförmig belastet ist, scheinen die Möglichkeiten der Anwendung dieses 
direkten Lösungsweges erschöpft, der einen geschlossenen Ausdruck für den statischen Anteil 
der Durchbiegung der Membranfläche ergibt. Bei Kreismembranen mit belasteten Teilgebieten, 
die von beliebigen Geradenzügen eingeschlossen werden, lassen sich keine einfachen partikulären 
Integrale der inhomogenen Differentialgleichung finden, die von der ersten Ordnung differen- 
zierbare Lösungsanteile liefern. Man erhält auf die angegebene Weise bei solchen stückweise 
gleichförmigen Belastungen geschlossene Ausdrücke nur für Näherungslösungen. 


Der Vorteil der Verwendung dieser Näherungen liegt im Wegfall des Gibbsschen Phäno- 
mens bei Reihenansätzen. Da wir wiederholt solche Näherungen benutzen, um den statischen 
Anteil der Lösung an den Übergängen vom belasteten in den unbelasteten Flächanteil ausreichend 
genau zu beschreiben, soll kurz darauf eingegangen werden. Für den konzentrischen Kreis- 
sektor gibt es eine Näherung für die statische Durchbiegung in der Form 


en BR e 2 3 (@ —r war ?) 
"A Te) 2 ee a+tr ar a 
a—r! a —r2 y (a? 72) a T 
ne er a er re ee re ?) 


es a —r2 v0 
+ a? sin 2 (9—9,) en In(?+@—2arcos(y—p))+ 5, °08 Ca) UBER 


BD,’ 
A 


2 2,2 
— 20 c0o®(@—9,)-+ a2cos2 (— 9,) o .. ’ 


+ 


mit den Festsetzungen 


0 für —-n <po <— m 0 für zn <o<0 
Balz BE eHitirt (<a, Gel 
0 für Drnen a re ee; 


Dieser stetige geschlossene Ausdruck ist keine von der ersten Ordnung stetig differenzier- 


bare Funktion, geht aber für 9, > ” in eine in den Stetigkeitsgebieten der Strömungsfunktion 


2 


hinreichend oft differenzierbare Funktion über. Für 9% -5 ist diese Funktion exakt der 


statische Teil der Lösung für eine auf einem Halbkreis gleichförmig belastete Kreismembran. 
Der zu dieser Näherung gehörige dynamische Durchbiegungsanteil gleicht formal dem dynami- 
schen Lösungsanteil für die Segmentlast. Die Reihenkoeffizienten dafür lauten wie folgt 


2 4 A — aa, Art n3 | 
Re z {—1+Jo(k,, Od} sinv9. + 22 Jar(k,, a) On sinvp + On sin 2790| - 


Diese Koeffizienten gehen für 9, — z über in die Reihenkoeffizienten des dynamischen Lösungs- 


anteils für die über einen Halbkreis gleichförmig belastete Membran. 


Der Umstand, daß die für die konzentrische Sektorlast gewonnene Näherung in die Lösung 
beim Halbkreis übergeht, kann als Kriterium ihrer Brauchbarkeit gelten. 


Manuskripteingang: 18. 11. 1958 


Anschrift: Dozent Dr. Max Hieke, Halle (Saale), Friedemann-Bach-Platz 6 
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Stress Distributions within Solids of Revolution 
By Avtar Singh 


In diesem Artikel wird das Problem der azxial-symmetrischen Druckverteilung bei halb-unendlichen 
iransversal isotropischen, durch ein oder zwei Kegel begrenzten Körper behandelt und formal mit Hilfe der 
M ellin-Transformation gelöst. ö e 

Als Beispiel wurde das Problem der Druckverteilung bei halbunendlichen Körpern unter gleichförmiger, 
auf ein Kreisgebiet ziehender Normalbelastung in Einzelheiten ausgearbeitet. Die Ergebnisse stimmen 
in dem speziellen Fall von Isotropie mit den von früher her bekannten überein. 


In this paper the problem of axially symmetric stress distribution within semi-infinite transversely iso- 
tropic solids bounded by one or two cones has been considered and a formal solution is obtained with the help 
of Mellin Transform. As an example, the problem of stress distribution within semi-infinite solids under 
a uniform normal loading applied over a circular area, is worked out in detail and the results in the special 
case of isotropy, agree with those previously known. 


Dans cet article le problöme de la distribution de pression azial-symmetrique dans le cas de corps demi- 
infinis et transversalement isotropes et limites par un ou deux cönes est traite et r&solu & l’aide de la trans- 
formation de M ellin. 

Comme exemple le probleme de la distribution de pression aux corps semi-infinis sous une charge nor- 
male et uniforme appliquee a un champs circulaire est Elabore en detail. Dans le cas special d’isotropie 
les resultats sont conformes a ceux connus du passe. 


B Hacrommefi pa6oTe pPaccMaTpuBaeTcH BOIPOC OCEeCHMMETPHYecKoToO PacıpenereHun 
MaBJIeHUA B IIOJIYÖeCKOHEeYHbIX, MHATOHAJIBHO U30TPOUHEIX TEeJIaX, OTPAHHYEHHEBIX OAHHM HJIH 
IByMA KOHYyCaMmH. Japaya pellaerTca POoPMAAIbBHO TPM TIOMOIIH HpeoöpasoBaHuun MernmHa. 
B xayecrbe IIpHMepa BoNpoc pacnıpenenleHun MaBlleHuA B IHOJIYÖeCKOHEeYHBIX TEeJaxX IPH 
YCJIIOBMWH HOPManbHOÄ MH PABHOMepHoÜ HAaTPy3KH, MelictBylomMefi Ha KpyToByIo IIIOIHaNB, 
pa3spa6aTbIBaeTtca MeTanmbH0. B, 4acTHoM CiIy4ae H30TPONHOCTH IOJIy4eHHLbIe Pe3ayJIbTAaTbı 
COBIANAIT C H3BECTHBIMH PaHee. 


1. Introduetion 


The problem of axially symmetric stress distribution in semi-infinite isotropice bodies of 
revolution, bounded by one or two cones, has been considered by Morgan [1] under quite general 
axially symmetric external loading conditions. The same problem has been extended here to 
the case when the material is transversely isotropic, with the axis of elastic symmetry, at each 
point, parallel to the axis of revolution. 


er 


< 


The semi-infinite conical shell of linearly varying thickness 


In section 2, asin Elliot‘s paper [2], the stresses and the displacements in a transversely 
isotropic material have been expressed in terms of two stress functions which satisfy two second 
order partial differential equations. These stress functions are analogous to the stress functions 
employed to define stress systems in isotropic materials. 
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In section 3, use has been made of the Mellin Transform to obtain formal solutions of the 
governing partial differential equations for the stress functions. It is assumed that the external 
loadings are such that their Mellin Transforms exist. Tranter [3], also used Mellin Transform 
for the determination of the state of stress within semi-infinite two dimensional wedges. 

If & and ß («x > ß) are the semi-vertex angles of the two bounding cones, then the solids 
(fig. 1) under consideration may be descriptively classified as follows 
Da>n/2>Pß>0 circular plate of infinite radius and linearly varying thickness, 

(ID) n/?22>&>ß>0 semi-infinite conical shell, 

(ID) z/2 >«&> 0, $ = 0 semi-infinite solid cone, 

(IV) a=n/2, ß = 0 semi-infinite solid, 

(V) « > n/2, ß = 0 semi-infinite solid with conical depression, 
(V\Da=n,n/2 >P > 0 semi-infinite solid with conical exclusion. 
In sections 5 and 6 case IV has been considered in some detail. 


2. Derivation of the fundamental equations 


The deformation being symmetrical with respect to the axis of revolution, cylindrical 
polar coordinates (r, ©,z) are employed with the z-axis parallel to the axis of symmetry at 
each point of the transversely isotropic material and coinciding with the axis of revolution. 
In these circumstances, the components of the displacement vector will be independent of the 
angle ©, so that we can assume 


UINT,2); = 0, DELL Rn se. 
The components of strain tensor are then ([7], page 56) 
ou u ow 
en eo , m d2’ 
ae (2,2): 


e == — —— B = == 0 
T2 dz ör ’ (CZ; ro 


Substituting these into the stress-strain relations for transversely isotropie materials we get 


w 


en ou u ö 
FIT ua tu 2) tn» 


— ou u ow 
90 = (Cy — 2 06) dr + Cu r + Cıa 22’ 


where c,1, Cs» C33s Ca and C,g arethe usual elastic constants. The stress equations of equilibrium 
are ([7], page 90) 
ir 180 2 m—60_ 


or r 00 02 r E | 
ro . 1200 0z 2r® | 
q = ara Due 
ör r 00 ja 02 2 r = aa 
0T2z 109z , &z L f2 en 


or 7 r 80 02 F 
Substituting (2.3) in (2.4) we see that the second of the equations (2.4) is identically satisfied 
and the other two reduce to 


u 1 u eu w 
a & Fr = re tt > 
as, 


Au 1 ou &w 1 ow &w 
(Ca ei a ef > taular % = Te = 
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Introducing the stress function ©fr, z) by means of the equations 


u ER NER ME RT 
r 
where k is a constant, we see that (2.5) are satisfied if 


&® 
HFidB+[lat klaut 12) Bor =0, 


&® 
[Ca + GC) + ku] ViD + Kin aa = 0: her Be 
where e an 
rn ET 
A ör + r ör 


A non-zero solution of these equations exists only when they are identical i. e., when 


K(&a + Cs) + Ca K Ca — 42 (sa 2.8 
A A Tr he y) } 9 MER 
Cı (a + C13) + Ko ( 
Elimination of k between these equations gives the following equation in »* 
Cr Cart + [2 Ca + C13) Cs — Er Cal 9? + Ca3 Ca = 0 ie On (A 


If »? and »3 are the roots of this equation and if ®, and Ö, satisfy the equations 


o | 
(+) =. = Ta... 22 Do 


then equations (2.7) are also satisfied. 
The displacements and the stresses in terms of ®, and ®, are given by the expressions: 


ae E ma > 2 (B,+D)+% E = +. k, =| 

50 — E ZT ni (+ D)+% E FE + k, =] em, 
zz=el,n—r0) i a8. + (kg Ca — 93 C1s) 2 i 
Fly ß + kı) ae +(1+k) Fe r9=-6:7=-0 


where k, and k, are the values of k corresponding to the values », and », of ». 


The roots »} and »; of equation (2.9) may be real or conjugate complex depending upon 
the elastic constants Cy,, Cag, Ca and Cjs. They are real, for example, in the case of magnesium 
but are conjugate complex for zinc. Thus we see from equations (2.8) that k, and k, are either 
real or conjugate complex. From equations (2.11) we infer, that for the displacements and the 
stresses to be real, the functions ®, and ®, must be either real or conjugate complex. This as- 
sumption is consistent with the equations (2.10). When »} and »3 are real positive, we will take 
»v, and », to be positive and when »? and »? are conjugate complex we will assume the real part 
of », and », to be positive. 


3. Making the substitutions 


” 


r=0_5Siny,, 
—_ — 0,c0osy, (dü=1,2) 
y; 


3.0) 
cosy =L, | 
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equations (2.10) are transformed into 


man =0 TER NTGE NN 


a 2 #l2u.ha. 1-48 
9 9 a ol, 4 


Multiplying throughout by 0% +! and integrating, from 0 to ©, with respect to o; we get 
MN DSDS DPD e0n een (83), 
where dashes denote differentiation with respect to Z,, and 


Dip, L) = f o?=1.@,(o,£;) de, (Mellin transform) a RE 


We get, on integration by parts 


; ro; r f PIOR ii, 
»+1 ae. F p oe R 
fe er AT IE je do, do; p ®; 
0 0 
MN R R 1272.30), 
oo, | 0, e 20, ” 
en ED: m : do; = ©" Pat dgie = R 
[e 08, do; 8; ’ f 0; o& do; od; s fe 00; 8; Ei p 2 ) 
0 0 0] h 
where we assume 
nd, g Hl) 
\ 8.6). 
» i ( 
E; 0%; —0d, 
as 0; — x 
The solution of equations (3.3) is 
Dir) A PB ONE Teer ar. 


where n =p—.1, P,„(£;) is the Legendre function of first kind of order n and argument Z,, Q,(&;) 
is the Legendre function of second kind of order n and argument Z,, and A, B; (i=1,2) are 
four constants of integration. 


Assume that 


u=m+Du,, v=w, +4, 
r=imt+trn, 69=-69,+60,  TE=TaHtzz, Br 3.8): 
RZearz+rz 
where 
0» od oD 0) 
I MN Du hehr,. mon 
— 10 [a 
r = (cn Vi — 209: ner D, +6 han: - 8:9), 
& 120 RD 
r nd 2 Cgg 94 uk 


and similarly for other stress components. 
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In terms of 0, £; these quantities become 


1-9) @ 
rn = & + — 9) u Hure r - z. 
{ 1-5) 2 
2 (CO act er Dr in 2m tt ar 2 


ö 
2) +9 


1) & 
Bu= [tn 2 + am gg Hat, en 


2utt ) 02 1 2 
— (a2 —4;) = N. ö£, + {ea — 200 — a)& + (u + re = 
(3.10), 
7 
+3. — cn +20) + (n—4ce—3a)} 0? | 
— _ [ra —H%Gel|, ‚ap A 
ö | vi I do: at 0? 06: ar 0; 60; 0%; 
Ban 1-40], 
a A R 
ee 2 (1 + Ki) Cya _ Cya 0? PR: (1 In 4) & 
Ty, = (1 — = Jule: el do? o? 06 
‚un @® 45 1-30 2], 
4 a Dr 
where 
a (= EI N EEE 


In spherical polar coordinates (og, y, ©), the stress components are given by the transformation 
relations 


oe =Trsiny+Fzsin2y-+ ZZcos’y, 
Yvp=FrTcos®y—fFzZsin2y + 7Zsiny, 
CIE ent... (3.12) 


En 1 nn mn n ee ‘ 
0y = 5 (rr —zz)sin2y+Fzcos2y. 


Assume again that 


— 


= 6h+0%, YP=-Whtim, 99-069, +09,, 0P=0oY+o0H, 
where (3.13), 


65; =fFr;sin®y+FZ%sin2y +Z2%c08y, (i=1,2) 
and similarly for other stress components. 


Substituting (3.10) into (3.13) we get, using the notation 


2 
Ki Ca — Pi Ca 


_g., “A+k) 


vi v 


hr a 
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— RD, 
09 = Kan + a a) siny + dcr py+ bl — sin 2yp]l 


+ Kar + eu) AD sinty + (1 drop hl EN sin pl 
1 88, 
+FRr@ ce )&l —&) sin?y + 24;6,(1— 2?) cos®y + 5,(1— 2) (1 — A)! sin 2y] — 0 00, ö£, 
1 28, 
+ (fa; + er Cs + (Cu — a) i}sin®y + d;(1 — LP) cos?y — b;L;(1 — LF)V? sin alle d0, 


+ [en + 260 — 34) +3, — cu) sinty —3d;L; (1 — Li) cos?y 


1 08, 


— 51 — 39) (1 — Ei)? sin 2 25 5a 


— : 3 ed, 
Fr = Mn + ad Eos y + sin NL sin 2p] 


1 29, 
+ Ka + (m — a) (1 — LP) cos?y +4; (1 — LH)? sin?y + bL;(1 —L})®? sin 2y] — 0? 08 


od, 
0; 00: O6; 


1 09, 
00; 


+2 EA Dos y+ 24Min 2NA— Ne sin2ylt 


(3-15). 
+ Kg; + Ca — 209 + (in — a)Eircos®y + d;,(1— 5) sin? y + 5,6; (1— EN)? sin aele 


+ Ken + 269 — 3a) + 3 (q; — C1) &i}&; copy — 3 d;L; (1 — Ci) sin?y 


1 08, 


+aA—3HA— sin dp] ng 


ob, 
oo; 


> 


oY = BG + (a; — Cı — Ad) i ——— +55 (1 — Gi)? cos 2») 


R si en 1 89, 
” le NT RN ee bi (1 — EP cos 2 En 
od; 
0; 00; 0%; 


1 
+ [&: 1 —&) {a — u — di} sin 2y + (1 — 225) (1 — CH)? b; cos 2 y] — 


od, 
0; 00; 


sin 2 . 


+ | 2 le +5 + (lu JE A1—Ed}- — b,&(l1 — SU cos |, 


‚sin 2 
+ | +20 30) +3 wir na: 2 


e% od, 
0; 0%; 


— (1 —32) (1 — EWR cos 2 2 


Denoting the Mellin transform by a superimposed ‘bar’ the following relations are ob- 
tained 
{20 = Sega da, de Pyi= SetvmeorT' "| 
- . (8.16). 
le: 00) = fe O&,or-1do, le - Ey = Serowıer-' de, | 
0 0 
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Multiplying equations (3.15) by 07 *' and ÄuigBEnE with respect to o; from . to ©, using equa- 
tions (3.16), we get 
{50 = IM DH) + PW+ DA —IMD, 
ee ii) DE+rE+D MIND, 3.17. 
let 80] = 19-9 + m PM) D+rW+DM— IND, 
ft 9 = MD + MPN) + PER + DM — MD 
RM, m. MR. MR. 
are the coefficients of 
ed, 2 Im BDO, 1 09; 
ie ah ar 


respectively, in the jth. equation of equations (3.15) with j=1,2,3,4 and i=1,2. The 
dashes denote differentiation with respect to {,. 
Boundary conditions: The prescribed boundary conditions are 


where 


Yy = gı(e) at y= x (outer Surface), 
— 9,(o) at y= f (inner Surface), 


=: u} . (3.18). 
ey=sile) at ya, | 
=s() at y=P ) 
Now suppose that 
Gu(p) = SE! Imde) er =! de, | 
(m = 1,2) st) 
S(Pp) = fi + Sm(e) 0? =! do . 
where g,„(o) and s„(o) are Mellin transformable external normal and shear loadings. 
From equations (3.1) and the relations 
r=osiny, 
Z=DCOSph u num wen 
cosy=L[ 
we arrive at the transformation relations 
b: = S - 
ee 
. (83.21), 
gg, = = Year) | 
or the equivalent relations 
Kun Gr 
ı—gı1—) ee SE A 


e= at -Hl—) 
Making use of these relations, the following equation is obtained 


le Dy oP ! do [® op y0 0” de + | 0° - Y Ya o” ' do 
= 
- [u ga dv RA AN 
ne } j Be > 
* [ü Zw}? Pride fl —2lı—H}]) (8.23). 


0 
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On the surface y = x equations (3.19) become, with the help of (3.23) 


ER micha. Rue 
GA =U—- RA MD. + Ba? Tea 
similarly, 8.24), 


LEER. a 
sp) =-U—-AL— N? K-owih-. + —- AA}? Me-övlh-. 
and on the surface y = ß 


Be TE € 2 EN 
GP) =-U—- BHA—M} ? Kr-pylb- +1 —- BA—H}? K-Pyll=a, 
?2+2 DEE Ge 
Sp) =-1—- AA} ? LK - Pb +1 —- 2A—M} ? Lei - 09-5 
where 
4,=4A{2+ (1 — Ay? = (1 + vjtano)t”®, 2 =6058 9. 2 843.26); 
A; = A{A2 +9 (1 — My = (1 + vi tan? BJ)", A = cos ß u a2 


Using Legendre differential equation to express the second derivative of functions P,(£;) 

. and Q,(£;) in terms of these functions and their first derivatives and then substituting equations 

(3.7) in second and fourth of equations (3.17) and the resulting equations in equations (3.24) 
and (3.25) the boundary conditions can be conveniently expressed as: 


G(p) = [Lı(p, 4, A191 kı) Pu(A) + Le(p, A, A, 91 Kı) Pr(A)llı Ar 
+ [L,(p, A, A191, kı) Qu(Ar) + Le(p, A, Ay 91 Ki) Qn(Ay)]lı Bi 
+ [Lı(p, A, Ay, 92, kg) Pa(A2) + Le(P, A, Ag v2, Ka) Pu(Az)] lz A, 
+ [L,, A, Ay, v2, Kg) QulA2) + Le(p; A, Ay v2 Ko) QnA)]oB : - - . - (3.28) » 
Gx(p) = [Ly(p, A, Ay, v2 kı) Pa(Ar) + L(p, A, Ay v3, Kı) Pa(Ay)] m, A, 
+ [Lı(p, A, Ay v1 Kı) QnAy) + Lx(p, A, Ay 91 Kı) Qn(Ar)] mı B; 
+ [Lı(p, A, As, v2, k;) P.(A,) + Ly(p, A, Az, v2, k,) P„(As)] m; A, 
+ [Lı(p, A, Ay 92, ko) Qu(Az) + Lx(p, A, Ay 92, Ko) QnAd)] m, B, . . . . (3.29). 
Sı(p) = [Mı(p, A, 4,91, Kı) Pu(d1) + My(p, A, A, v2 kı) Pa(Ay)]ı Ar 
+ [Mı@p, A, A, v1, Kı) Qu(A1) + Map, A, Ay v1 Kı) Qn(A)lı Bı 
+ [M;(p, A, Aa, v2, k,) Pa(A2) + Ma(p, A, Ay v2, ka) Pr(A2)] I, A, 
+ [Mı(p, A, Ag, v2, ka) Qu(A2) + My(p, A, Ay v2 Kı) QnA)lEB, - . . . (8.30). 
S,(p) = [Mı(p, A, Ay v1 ky) Pu(Aı) + Map, A, Ay 91 Kı) Pr(Aı)] mı A, 
+ [Mı(p, A, Ay v1, Kı) Qn(Ar) + Ma(p, A, Ay v1, Kı) On(Ay)] mı B,' 
+ [M,#, 4, A, v2 ka) Pa(A2) + Ma(p; A, Ay v2 ka) Pu(A,)] m, A, 


+ [M,;(p, A, As v2 K,) Q(A,) + My(p, A, Ay v2 Ka) QAylmy B,-.: » » 831) 
„ where 
Dr Pr 
We Al dt leiten (3.32), 
and e 
_ 2 an ® 
Lp&wk)=; ahnt (fi5 p fi ’ 
(2) (2) Dr 1 (2) 
L;(p, &, ie) Zpip + fi ha —ı- a Ji2 
1—£& 
SE (3.33), 
26 


M,(p; 6: C Vi» k;) — 1 


[ 
+ No, | 
Map, L, Cs Yis k;) = p En 1) y a 2 Fr a el | 


and for instance 


Pa = | 7 >|, BR Pr Erkn, 
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The four equations (3.28—3.31) determine the four constants A,, A, B, B, and thus the 
Mellin transform ©; of ©, is completely known. Once ®, is known, the stress function ®, can 
be determined from the Mellin’s inversion formula 


e+i® 


(0,6) = or -Din.Eran een Di . (8.35). 


Zei ni 

The imaginary number i = y— 1 occuring in this SUB UOHE has no chance of being con- 
fused with the subscript {. 

However, more interest lies in the determination of stresses themselves, a further appli- 
cation of Mellin’s inversion formula gives 


I i 


Da > 
RT Eh vor 0; ?{ot - 004 dp, 
ri 

1 N 
ee ee —p 2. ee 
7 72 Fe | 0: Piel pp dp 
ee str a 
eg 1 s een 
d:00=-,_; 0, ?\o?- 00, dp, 
oo; | retten) dp 
} ee Se : > 


where the barred quantities in equations (3.36) are given by equations (3.17). This gives, for- 
mally, the complete solution of the problem. 


4. Stress distributions within semi-infinite transversely isotropie solid cones 


In the case of a solid cone (3 = 0) only two boundary conditions can be prescribed namely 


yy=y oe andy = sı(o) on the outer surface y = a. Each of the functions ®, and ©, contains 
two constants of integration, two out of these four constants are, therefore, to be disposed of. 
This can be achieved by applying the condition that the stresses must be finite on the z-axis, 
where; ={=1(y=0). This requires that the constants B, = B, = as the function Q,(d,) 
and its derivative with respect to £,; has singularity at £; = 1. On the other hand the function 
P,(&) 7» sis derivative with respect to £; is regular on the z-axis. The functions ©, are then 
given by 

Dip, &,) — A; P.(&) a ee TE (4.1). 


. Stress-distributions within semi-inflinite {iransversely isotropie solids subjeeted to Mellin 
transformable axially symmetrie loadings 


In this case the angle & is taken to be z/2 and the functions ®, are given by (4.1). 


In this case it is convenient to work with equations (3.10). On taking the Mellin transform 
of 05-22, and 0: r2, we have 
{22} = 41 MH SA HAap+I)d+rP +9 — 18], 
Hr Fl HIER-WLY—-AP+ ING — LI — Hi] 


Equations (5.1) when evaluated at the boundary of the semiinfinite solid, where © =&; —0 
(y=aj2) give 


(5.1). 


{0° 22.0,D} = — dp PDA, un u 
{0 7240, Pd} = —bi(p +1) Pi(0)4, 2 12 Re 
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The values of the Legendre function P,„(£;) and its first derivative at £; = O0 can be expressed 
as ([5], page 144) 


Vr — 2yr 
P Fr — 2 — er: — 
(0) r 2+n r{! Ren ’ P,„(0) r 1+ n 7 n . (5.4). 
2 we ae 
Substituting in (5.2) and (5.3) and remembering that n = p — 1, we get 
ee d; p? m 
fo? . 22,(0, p)}} = I N .2(3:5); 
Te 
2 
— 2b; 1) yr 
ee LEE matt RR (5.6) 


Substituting (5.5) and (56) in equations (3.24) with A, = A, = 0, the following boundary con- 
ditions are obtained 


: pr Yr 
En A EEE EEE 
IE pP 2—p 
jaren 
A RE Anh  Asbe)ie ehren DR 
ee 


When solved for A, and A,, these two equations give 


a Hp) GP) bs — BP) SD) de, 


Di eb.d, 
33) 
1, HP) GP) dı — Alp) Sı(p) dı 
i bi d,— b,.d, 
where ee > . 
p BR p —=E 
Kae) „,_rtzN 
H(p) = - — ee, HA, = u IE). 
Now we find Z2 at the ee axis, wreery=0,l=-6 = 
using the result Ps) = en 
Making use of equations (3.23) and (5.11) with £, = 1, we obtain Per 
{o® - z2(1, p)} =pp+Y)mrdA +nNt’d, As} 
1 
RL ar Tr [dv (HG b, — MH, S,:d,) 
WENN, 5 de HE STR ol 


The stress ZZ on the axis is obtained by applying the Mellin inversion formula to (5.12) 


| BP ZELL DIN aha anne we 10-18): 
6. As an example we work out the case for a semi-infinite solid subjeeted to a uniform normal 
loading P per unit area, applied over a eireular area of radius R 
The boundary conditions now become 
j P, eur‘, 
\0 : R= 0-60 
or Da NE ea EEE NT 


Ile) = 
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Equations (3.19) then give 
p+2 
GP) = ——ı S(y=0: 2.0.0. 3 ae 


Substitution in equation (5.12) and the resulting equation in (5.13) gives 


ce + ioc 
a 2 (b,d,y +? — da, #}) et © Nie 
I ee | | r{fP&®\rt ZP\eir) ap. 
R) 2 Aniyn (b, d, —b, d,) | p+2 2 2 


(6.4) 


CGauchy’s theorem of residues can be made use of to evaluate this line integral, taken in the 
complex p-plane. The Gamma function /‘(p) can be represented thus ([4], page 207) 


. . oo (— I 

KXp)=# A :- 
P=Fo+ 2 arm ) 
where F(p) is an integral function. The poles are therefore givenbyp= —n(n=0(,1,2,...). 
Therefore the singularities of /'(— p/2) and I’ ers aregivenbyp=2nandp=—(2n-+3) 
withn=0,1,2,.... Thus the integrand in (6.4) has simple poles at p = — 2, — (@ n +3), 
2n (n=0,1,2,.... For —2<c<0Oand (o/R)< 1, and taking a contour as in reference 

[1], the sum of the residues of the integrand /atp=—2 and —(2n+35) gives 

c+ix 

ag Idp = Ya (e/RP [(.dı did) + did, (e/Rv){l + (e/Rv | (6.6), 


e—io 


— bad, (e[Rv){l + (e/Rv)y"P], o/R<1 
The residue at the poles p = 2 n is zero. The same result (6.6) is obtained for o[/R>1. 
Substituting (6.6) in (6.4) we get 


1 Ib, d, o/R b.d, o/R 


b,dh, —bid, | va u ba 62 02 of... (6.7). 
1 me Ss 
a al 


In the case of isotropic material, =, =l,,=k=1,d4d, =4d, = (a — „and b, = b, = 2cu- 
The equation (6.4) in this limiting case reduces to 


zz(1,0=P|i1+ 


200 


a (p +3) nr x rl - e) 
(o/R)* - ZZ(1, 0) = ve — (o/R)-rdp . . . (6.8), 
r 2rniyn (p + 2) (iR) p NR, 
giving > 
En h 0° 
z2(1,0) = “are Ep am) TEN 


‚ Ihe same result is, of course, obtained from equation (6.7). This result agrees with that 
obtained by Timoshenko [6] except for a change in sign of P due to the difference in the direc- 
tion of the applied load. 

Stress distributions are thus obtained under quite general conditions of external loading. 
The displacements are determined making use of the stress-strain relations (2.3). 


’ The analysis developed here might find an application in some problems of technical 
ımportance such as the determination of the state of stressin conical pivots and in conical domes. 


The author wishes to thank Dr. D. N. Mitra for guidance and helpful suggestions during 


the BERND of this paper, and to Prof. Dr. B.R. Seth for his encouragement during the 
work. 
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KLEINE MITTEILUNGEN 


Zur Berechnung der reibungslosen inkompres- die trigonometrischen Veränderlichen g bzw. y ein- 
siblen Strömung für ebene Schaufelgitter geführt: 


Einleitung Egal (1 — cos p), 


Die Schaufelgitter bilden eines der grundlegenden ze R| (1). 
Elemente der Strömungsmaschinen, und auf Grund Fl 
dieses Umstandes arbeitet man ununterbrochen an ä a 
verschiedenen Methoden der Berechnung von aero- 
dynamischen Eigenschaften der Schaufelgitter. Für Die Geschwindigkeit, die an der Sehne der Grund- 
die Berechnung von nicht allzu dicken und mäßig schaufel von den an der (— v)-ten und »-ten Schaufel 
gewölbten Schaufelgittern ist die 
Singularitätenmethode, bei der 
die Schaufeln durch eine konti- 
nuierliche Verteilung von Wirbeln 
und Quellen ersetzt werden, am 
meisten geeignet. Diese Methode 
wurde von mehreren Forschern 
ausgearbeitet und vor allem in den 
Arbeiten von H. Schlichting [1] 
und N. Scholz [2] einer leichten 
und raschen numerischen Berech- 
nung angepaßt, freilich unter der 
Bedingung, daß die Koeffizienten 
für die Berechnung induzierter Ge- 
schwindigkeit, sog. Abwindfunk- 
tionen, bekannt sind. Deshalb hat 
auch H. Schlichting seiner Ar- 
beit [1] noch die Tafeln dieser 
Funktionen und zwar für Teilungs- 
verhältnisse 2/1 = 0,5; 0,75; 1,0; 
1,25; 1,5 und 2,0 sowie Staffelungs- 
winkel A = 0°, 15°, 30°, 45°, 60° 
und 75° beigeschlossen, und da- 
mit hat er eigentlich eine breite 
Anwendung der Singularitäten- 
methode zur Berechnung von 
Schaufelgittern ermöglicht. Aber in der Praxis be- verteilten Wirbeln induziert wird, ist mit dem Aus- 
gegnet man oft den Schaufelgittern mit einem anderen druck: 
Teilungsverhältnis und Staffelungswinkel als den von waslE) = Us — it 
Schlichting tabellarisierten, und deshalb wird an / Re 
dieser Stelle eine Methode für eine einfache Berech- i [ 
2n 


1 
= (1l—cosy), 


nung dieser Koeffizienten für beliebige Teilungsver- 
hältnisse und Staffelungswinkel geliefert. 


Kr 

un) 
l 

> 


« 


Kim 
r 


Zur Ableitung hat man dieselbe Bezeichnung wie 
H. Schlichting [1] verwendet. Gleichfalls werden 


gr 


anstatt der dimensionslosen Koordinaten E bzw. & 


r Rn l R en 
Theorie er ern = y(E) dE 
2 y — ei 


seben. 


nu 


Ie 


) 


le 
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In den Ausdruck (2) setzt man die trigonometri- 
schen Veränderlichen 9, y ein, drückt die Zirkula- 
tionsdichte y(&) durch die trigonometrische Reihe: 
y(& =2U_|Asctgy +Arainy + Arein2yt ... 


(3) 
aus und drückt das Teilungsverhältnis und den 
Staffelungswinkel durch einen einzigen komplexen 


Parameter x =i en eriAaus. Nach unbedeutenden 


Umformungen ergibt sich: 


Te 
Was Han na N TE 
U% Us U n cosp —cosy va 
ü) 
E i 
cos p — cosy— va 


R ja + c08 %p) lc (1— eos 2%) 


+ Ar eo na) dr N) 


Gleichfalls wird noch eine Bezeichnung für die Inte- 


grale 
nr 


1 cosny 
) mr - | 
Dn(p) en re 

ö 

5 
nl 1 | DER Oi 
In cos p — cosy —vı 
0) 


eingeführt. 
Jetzt kann man die Beiträge der (—v)-ten und v-ten 
Schaufel zur Abwindfunktion ausdrücken: 


Aylfyn + %9yn) = Da-ı(P) — Dar(P) ; 
EEE PR (6), 
4, tim) =2Die)+2De@) ----..M. 


Es bleibt noch übrig, die beiden Integrale in der 
Gl. (5) zu bestimmen. Mit Bezug darauf, daß der 
Kosinus eine gerade Funktion ist, kann man schreiben: 


Pr 1 
1 cos n 
Date) = zu | —_ nr dy 
r 


c0o8 pP — cosy va 


Bi 
un 2 cosny j g 
— —— Me RlE)* 
RE TTTE i. (8) 
er 
Wird von uns statt der Veränderlichen y die kom- 
plexe Veränderliche w eingeführt: 
Du... EB, Dr 
so gilt 
co3np = ex un -} 4 dy . 10 
— 3 | rk ' ww (19), 


und die Integrale (8) gehen in die Integrale 


l 
u ER N Fr 1 
BR Er we —2(cosp +vo)w-I ® 


e|=1 


> wn — — 
ı wn 


r 4 w? — 2 (cos Po—v o)w-+ 1 in) 
lv) = 
3 a h - un dw 
2n w — 2 (cosp +vo)w+1 
w|=1 
ı URAN un dw 1 
an we —2 (osp—vo)w-+1 En 


über, wobei der Integrationsweg der im positiven 
Sinn durchgelaufene Einheitskreis ist. Die Wurzeln 


des Nenners im Integranden des ersten Integrals sind 
3 = cosp+vaFYra?+2vacosp—sin?p (12), 


- wobei 


e|>1, w|>1, 


und die Wurzeln des Nenners im Integranden des 
zweiten Integrals sind 


wa 0089 —v2+ (v! 22—2 va cos p—sin?@ (13), 
wobei 


| <ı, jJol>1. 


Innerhalb des Integrationsweges haben also die 
beiden Integranden je einen Pol der ersten Ordnung, 
so daß im Sinne des Residuumsatzes gilt: 
| u uw 

w—W W— Us 
Werden diese Ausdrücke in die Gl. (6) eingesetzt, 
ergibt sich 


27 (fyn "Fr i gyn) 


... (14). 


Day) = — 


l 


Setzt man in die Gleichung (7) ein, ergibt sich 


Sue 


WW — Wg 


9 1+% 


Ws — uy 


Ayltyo + I gy3) =—2 


oder 


3 l1+w 
A, (fyo +ty) = ! ——— 
\r?a2+2»xcosp — sin? 
£ 1+w 
ra — 2va&cosp — sin?p 
(16). 
Für ziemlich große » kann man die Ausdrücke 
für w, und w, in Potenzreihen nach den Potenzen 


1 f 
— entwickeln: 
va 


j > sm n2 17 
_ = CcOoSso sın“ o\/ 
Wa cos ra Fva (1 En zu q ein ) Ei 


1 
Yale ee 
RN va 2yl a? 


sinp 1 (5 9 __ sin? 3 
5) 


I ( p 2) 5 [22 sin? pt 
le nn | | 


va 22a? sIyx Iy2 a8 


_ 7 /cosp __ sin? 9\5 
gl 7 | + . (17), 


vX  Iy2 a2 


wo die oberen Zeichen für w, und die unteren für Ws 
, 2 r 

Geltung haben. Nach einer einfachen Umformung 
bekommt man 
Alhı +8 Iyı) wı 7 u, = coS ® 


Wr 


l 3 1 \ 
4 7 cosPp + FR cos 3 o) 


ST l 
Erd 2 a 17 


1 /245 nr 7 
45 (35 0089 an 150839 . 15 °0s5P) 
Yet RN 1323 315 

„10 | Fra Hr 32 c08 3 ’) (18). 
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Der gesuchte Funktionswert (f,, + # 9,,) ist dann 
durch die Summe der Ausdrücke (18) für » von 1 
bis © gegeben: 


ö = i 1} 
In + im = PY Iıllyı +39) = 362) cos p 
v—1 


di 3 1 
ze (20059 +7 0839) 


1 25 15 1 
= DE 


l 245 49 7 
— 8: (32 0° P + 75008 39 7 e 


1 1323 315 
— m &10) : (sr cos p +30) ae, 
(19) 
wo &(n) die Riemannschen Zetafunktionen 
4 1 
Sn) = en 


sind. 

Die Reihe (19) konvergiert für |x| = 2-1/l > 2. Um 
die Entwicklung (19) auch auf die dichteren Schaufel- 
gitter anwenden zu können und dabei noch eine gute 
Konvergenz gewährleistet zu haben, wird man den 
Einfluß einiger naher Schaufeln nach den genauen 
Formeln (15) berechnen und erst für die restlichen 
die Entwicklung (19) verwenden. Wie später erwähnt 
wird, genügt es, wenn man nach der exakten For- 
mel (15) nur den Einfluß von einem oder zwei be- 
nachbarten Schaufelpaaren berechnet und nur für 
sehr dichte Schaufelgitter ist es erforderlich, nach den 
exakten Formeln den Einfluß von drei Schaufelpaaren 
zu berechnen, um die ausreichend genauen Werte zu 
erzielen. Die entsprechenden Entwicklungen ergeben 
sich aus der Entwicklung (19), falls man darin anstatt 


von &(n) fortschreitend £,(n), &3(n) bzw. £,(n) schreibt, 
wo 


. (20) 
ist. 
Ganz ähnlich geht man auch bei der Berechnung 


restlicher Abwindfunktionen vor. Darum werden im 
weiteren nur die endlichen Formeln angeführt. 


I Ayla +8 gy0) = z (0,6449 — 1,2899 cos p) 


+ en (0,1852 — 0,3705 cos@ + 0,1235 cos 2 
& 
— 0,0412 cos 3 p) 
1 
E= =. (0,108 — 0,217 cos 9 + 0,108 cos 2 @ 


— 0,054 cos3@ + 0,011 cos4 
— 0,002 cos 5 p) 


un m (0,078 — 0,156 cos + 0,094 cos2@ 
— 0,056 cos3 9 + 0,019 cos 4 @) 
+ > (0,06 — 0,12 cos @ + 0,08 cos 2 @ 
— 0,05 cos3@ + 0,02 cos49p) ; (21), 
34, (hi + tmı) = — = - 0,6449 cos 
= (0,1235 cos 9 + 0,0206 cos 3 p) 


— . (0,054 cos + 0,016 cos 3 

+ 0,001 cos 5 ) 

— <= (0,031 cos + 0,012 c0s 3 p) 
& 


— m (0,02.c0s @ + 0,01 c0s3 p) ; (22), 
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= N 
Z Aula + ip) = + 25 0,3225 
+ z (0,0823 + 0,0617 cos 2) 
+ = (0,041 + 0,043 cos 29 + 0,005 cos 4 p) 
[0 4 


E= G (0,025 + 0,031 cos 2 9) 


+ (0,02 + 0,02 0029); . (23), 
34 ee 0,0617 c08 p 

2 (0,049 008 P + OHM E65 89) 

u =: (0,037 cos@ + 0,013 cos 3 9) 

— — (0,03 cos + 0,01 cos 3 p) » (24) 


Die Formeln für die Summen mit den ausgelassenen 
zwei oder drei benachbarten Schaufelpaaren ergeben 
sich aus den Formeln (21)—(24), wenn man die Aus- 
&s(n) bzw: A) 
&a(n) "&(n) 
multipliziert. Die Werte dieser Koeffizienten werden 
in der Tabelle 1 angeführt. 


drücke bei Es mit den Koeffizienten 
& 


Tabelle 1 
n.. 1.2, 
Gn)/&n) | 0,6124 | 0,2408 0,0991 | 0,0420 0,0181 
tum)/&n) | 0,4401 | 0,0908 | 0,0200 | 0,0046 | 0,0011 


Ergebnisse 


Zur Erleichterung der numerischen Berechnung 
werden noch die Formeln für die dem sog. Dreiviertel- 
sehnen-Theorem entsprechenden drei Punkte an- 


4 Efar BER: 
geführt, d.i. für x/l = B’I’R 
Für 
3 
zjl = 75 ee! 


berechnet man 4, (fyn + i 9yn) nach den Formeln (15) 
und (16) und dazu addiert man den Einfluß der rest- 
lichen Schaufelpaare: 


- — 0,0206 tt cos 4A + 0,007 7° cos 64 


ılM8 
tv BR 
“ 
en 
> 
— 
|) 


— 0,002 78 cos8A; 
3 A, fyı = 0,3225 7? cos 24 — 0,0412 1? cos4 A 
2 


-+ 0,011 76 cos 6A — 0,004 73 cos8A ; 


I Ay hype = — 0,3225 7? cos 2A + 0,0515 7! cos 4A 
he — 0,016 76 cos 6% + 0,006 73 cos8A; 
> Ay fys = — 0,0309 1 cos 44 + 0,014 16 cos6A 


— 0,006 73 cos 8A . . (25). 


er 3 
wo noch die Bezeichnung =, ; eingeführt wurde; 


t 
l 


Dr 


für’ > > 0,65 
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berechnet man nach den genauen Formeln j 7 
3 == — 
5 a 12 
A, (hy +19)» At 9) 
3 tt zı 
und für die Berechnung des Einflusses aller restlichen Er 
Schaufeln gilt: 3 Ayfn = — 0,8599 12 cos 24 + 0,1105 ri cos 4A 
=? 
3 Ay fu = — 0,0050 7 cos 4A + 0,0007 7° cos 64 — 0,026 7° cos 6% + 0,007 7% 00584; 
Fön —_ 0,0001 73 cos 8A ; ZA, fr, = — 0,1075 1? cos 24 + 0,0107 rt cos&A_ 
’ ’ a 2 
x — 0,002 19 cos 64 + 0,001 7 cos 84; 
3 4, fyı = 0,1975 7? cos 24 — 0,0099 7? cos 41 4 
RE YA, fy = — 0,3225 7? cos 24 + 0,0240 rt cos 44 
— 0,0011 7% cos 6 A — 0,0002 73 cos 8%; ‚=? 
— 0,004 78 cos 64 + 0,001 7? cos8A; 
{ =. - ? ri cos 4} © 
EEE REDE PEHRUNIERF IE 3 Ay fin = 0,0103 74 cos 4A — 0,003 7° cos 64 
a 6 8 f R v=2 
0,0016 76 cos 6A + 0,0003 73 cos 8%; + 0,001 78 008 84; 
3 24 fys = — 0,0074 ri cos 4% + 0,0013 7 cos 67 L > tl = 0,65: 
= 
ef. s . (6) 
DENE 2) in = — 0,5266 12 008 24 + 0,0266 rt 000 4A | 
für 0,65 >t/l > 0,46 | 
0,0026 7? cos 64 + 0,0003 7? cos8% ; 
muß man nach den genauen Formeln den Einfluß von r 
den ersten drei Schaufelpaaren berechnen, indem für 57 Am 0,0658 tr? cos 2% + 0,0026 r!i cos4 4 
alle restlichen die Gleichungen gelten: dt 
- 0,0002 7° cos6% ; 
= in rs , S 
ae a a er en 8 3° 3 4A, I, = — 0,1975 t? cos 2% + 0,0058 rt cos 4A 
„=3 


0,0002 7 j} 
—(, j6°088 & 


2 4 
Ay fr = 0,2838 n cos 2A — 0,0150 T cos 4} 


ıMe& 
En 


v 


6 8 
+ 0,0018 z cos 6 A — 0,0003 2 cos 84 


10 
+ 0,0001 35 cos 104: 


= 12 : a L 
F A, ty, = — 0,2838 ; cos 2A + 0,0187 y hi 4) 
= 0) 0026 5, — 086% +4, 0005 7 „cos si 
—(, 0001 ;, , cos 104; 
zi 76 
H fy, = — 01,112 1 cos 4) + 0,0022 r cos 6A 
- 0,0004 nn cos 84 + 0,0001 m cos 104 


(27). 


Die Formeln für & 4, 9y„ ergeben sich aus den ent- 
sprechenden Formeln für 2A, f,, wenn man darin 
sin mA anstatt von cos mA schreibt. Berechnet man 
die Werte der Abwindfunktionen nach diesen For- 
meln, ist der resultierende Fehler kleiner als 210-4, 
Wenn nicht so große Genauigkeit der Berechnung 
gefordert wird, kann man in den erwähnten Ent- 
wicklungen eine kleinere Anzahl der Glieder nehmen 
und auch den Anwendungsbereich der Formeln in der 
Richtung kleinerer Teilungsverhältnisse verschieben. 


e 7 

Für »/l = BE und z/l = 
und deshalb werden nur die Formeln für die Sum- 
men erwähnt. 


1l 
12 gilt dasselbe Verfahren 


— 0,0004 1% cos 6A : 


3 A, fy, = 0,0025 rt? cos 4 4 — 0,0003 7° cos 64; 
=3 


‚65 > t/l > 0,46: 
} T? tt 
2 4A, ty = — 0,7569 5 cos 24 + 0,0402 ‘ cos4/ 
7 - 
7$ 
- 0,0042 z cos 64 +0, 0005 | - cos By 
z!0 > 
- 0,0001 33 cos 102: 
x T? r! 
I 2 ds Iı - 0,0946 = cos 2% + 0,0039 z cos 4, 
gm 3 
” L 
0,0003 g 608 643 
cin . ig 
I, /y, — 0,2838 , cos 2% + 0,0087 g 008 47 
- 0,0006 z cos 64 + 0,0001 es cos 8A; 
u r! . 1® 
& 4A, ty, = 0,0037 7Y cos 44 — 0,0005 5 c08 64 
+0, 0001 16 ni 3A, 
LE 
B ee a 
tl | 
Bi 1,7198 7? cos 2A + 0,5343 rt! cos 4% 


— 0,311 7% cos 6A + 0,214 73 cos 8% 
— 0,16 r!0 cos 104 ; 


ZAMM 39 (1959) Heft 12 i ; Kleine Mitteilungen 499 
© 0,65 > t/I > 0,46: 
3 A,fr, = — 0,5374 72 cos 2A + 0,0991 tt cos 4A ent 
42 © T> T# R 
— 0,041 78 cos 6% + 0,022 73 cos 8% er Ay fa = — 1,5137 5,008 2A + 0,1941 go 4 
Ps n 8 
— 0,01 r10 cos 10%; — 0,0498 = cos 6% + 0,0157 In cos 8A 
x. 
3 A,fy = — 0,3225 7? 2. ‚106 - 10 
zit T? cos2A + 0,1063 rt cos 4A — 0,0055 55. cos 104: 
 — 0,054 76 cos 6% + 0,032 73 cos 8% i > 
[} T 4 
0,02 120 cos 104: „2, A, = — 0,4730 > cos2i + 0,0360 7 cos4/ 
x 76 a 13 
I Ay hy, = 0,0515 Ti cos 42 — 0,039 7% cos 6A = a 
+ 0,027 73 cos 8A — 0,02 r!0 cos 104. 0,0005 35 cos 102: 
1 > 0,65: x ? 
NE =9,68: 3 Art = 0,2838 Z 00824 + 0,0386 e cos 4A 
x v—4 R 
5 A, fr = — 1,0532 T? cos 2A + 0,1287 tt cos 4} 16 13 
v—3 — 0,0086 2 08 64 + 0,0024 Tg cos 81 
— 0,0309 7% cos 6% + 0,0090 73 cos 8% 
_ ro 
— 0,0029 r10 cos 104; EEE 
n3 oo m Tr! Br ©) a aa 
3 A, fu = — 0,3291 72 008 2A + 0,0239 72 cos 4 ee N 
3 zZ 
— 0,0041 7° cos 6% — 0,0009 7? cos 84 + 0,0020 © cos 81 
— 0,0002 710 cos 104 ; 710 
— 0,0007 32 cos 10%. 
3 A, fy, = — 0,1975 72 cos 2% + 0,0256 7? cos 4A ; 
v_ Beispiele 
— 0,0053 7% cos 64 + 0,0013 73 cos 8% Mit erwähnter Methode berechnet man die Werte 
der Abwindfunktionen für die Schaufelgitter mit den 
_ 4 110 cos 104 ; 2 : =. ö 
ee Teilungsverhältnissen t/! = /2 und {/l = und mit 
B> Ay fr = 0,0124 7! cos 4A — 0,0038 1° cos 64 dem Staffelungswinkel von 4 = 30°. Die Berech- 
‚»=3 


+ 0,0012 73 cos 8A — 0,0004 r!P cos 104 - 


nungsergebnisse für t/!— 2 sind in den Tabellen 
2—5 zusammengestellt: 


Tabelle 2 
x/l w, Ws Ar fyı | Dial Iyı fyı A, 9yı | 34, 91 Iyı 
| IE 2 
B 0,0952 0,1282 | —0,0612+i0,1757 | 0,0340. 0,0205] 0,0545| 0,0475 0,0344] 0,0819 
15 | 007691 0,1609 | —0,088141 0,1446, —0,0112 | —0,0068 | —0,0180 | —0,0163 | —0,0115 | —0,0278 
| | | | 
. 0,0413 0,1875 | —0,0989-Hi 0,1128 | —0,0576 | —-0,0343 | —0,0919 | —0,0747 | —0,0569 | — 0,1316 
Tabelle 3 
er? | 14 Labs - 
z/1|- a I  — —_ I An I|Shm| mo | Am |EZAm| m 
Vx?+2xcosp — sin? IVa®— 2x cosp— sin? | 2 | 2 
2| 0,1669 0,3071 | 0,1657—: 0,3038 0,0013 0,0002 | 0,0015 | —-0,0033 | —0,0003 | —0,0036 
12 | | | 
5 0,1027—1 0,3665 | 0,1938 0,2400 | —0,0911 | —-0,0546 | —0,1457 | —-0,1265 | —0,0916 | —0,2181 
| - 0,0093 0,3930 | 0,1990 0,1848 | —0,1897 | 0,1111 | —0,3008 | —0,2082 | —0,1790 | —0,3872 
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Tabelle 4 
mm Te 
z/l uw | uw; Ay fr | 24, Ir | Irı | Aı 99, z4r Iyı Iyı 


een, ze nissen kt — 


fl j 
—0,0074—i 0,0244 
—0,0200—i 0,0247 | 


| 


—0,0335—i 0,0155 


SIE Bl wie 


Tabelle 5 


0,0271: 0,0215 | —0,0345 | —0,0205 | —0,0550 | —0,0459 | —0,0342 | —0,0801 
j 
—0,0131—i 0,0255 | —0,0331  —0,0203 | —0,0535 | —0,0502 | —0,0346 | —0,0848 


' 0,0029—i 0,0223 | —0,0364 | —0,0209 | —0,0573 | —0,0378 —0,0335 | —0,0713 


D) A,fı I A, 91 240m I9y 
z 0,0038 —i 0,0014 | 0,0054—i 0,0034 0,0016 | 0,0002 0,0018 | —0,0048 | —0,0004 | —0,0052 
. —0,0055 +: 0,0013 | 0,0048 +1 0,0004 | —0,0007 | —0,0001 | —0,0007 0,0017 | 0,0001 0,0018 
= —0,0043-+i 0,0056 | 0,0028-+i 0,0019  —0,0015 | —0,0003 | —0,0018 0,0075 | 0,0007 0,0082 | 


Bei dem zweiten Schaufelgitter, i/I = m ‚A305 


muß man den Einfluß von zwei Schaufelpaaren nach 
den exakten Formeln berechnen, es ist aber dabei zu 
beachten, daß das zweite Schaufelpaar bei diesem 


Schaufelgitter zu den Abwindfunktionen denselben 
Beitrag gibt wie das erste Schaufelpaar bei dem 
vorangehenden Schaufelgitter (t/l = [2 ). Die Ergeb- 
nisse der Berechnungen sind in den Tabellen 6—9 
zusammengestellt. 


Tabelle 
Fa er ZAm| I EraRr ZAm| m 
= 0,0030 | 0,0013 0,0006 0,0048 —0,0414 | —0,0033 | —0,0010 | —0,0457 
= —0,3749 | —0,0911 | —0,1347 | —0,6007 —0,3357 | —0,1265 | —0,2223 | —0,6845 
z —0,6151 | —0,1897 | —0,2759 | —1,0806 | —0,3609 —0,2082 | —0,4286 | —0,9977 
Tabelle 7 
xl Aıfyı | An | Z Art | fr | dı 9yı A, 9yı 34, 97, Iyı 
R 
n 0,1381 | 0,0340 0,0505 0,2226 | 0,1304 0,0475 | 0,0834 0,2613 
> 0,0474 | —0,0112 | —0,0168 | —0,0754 —0,0476 | —0,0163 | —0,0280 | —0,0918 
ee —0,2181 | —0,0576 | 0,0848 0,3605 | —0,1823 | —0,0747 | —0,1374 | —0,3944 
Tabelle 8 
z/l A, fr | Ay fys | 3 Ayfy | fr | A, 9, A, gr; | 34, Iys gys 


| 


| 
a ep er 
| | | 


r —0,1358 | —0,0345 | —0,0507 | 0,2210 0,1137 | —0,0459 , —0,0829 | 0,2425 
7 i 

1“ — 0,1404 | —0,0831 | —0,0501 | —0,2235 | —0,1598 | —0,0502 | —0,0843 | —0,2943 
11 


—0,1184 | —0,0364 | —0,0520 | —0,2067 


—0,0453 | —0,0378 | —0,0800 | —0,1631 
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A | A At | ZA | Ban | Wr ZAm| m 
IT  ———— jpc Ol 
> —0,0037 | 0,0016 0,0008 | —0,0013 —0,0527 | —0,0048 | —0,0016 | —0,0591 
5 —0,0015 | —0,0007 | —0,0003 Re, 0,0215 0,0017 0,0005 0,0237 
| = 0,0207 | —0,0015 | —0,0011 0,0181 0,0572 0,0075 0,0026 0,0673 
en [2] N. Scholz, Strömungsuntersuchungen an Schaufelgittern. Teil. 


[1] H. Schlichting, Berechnung der reibungslosen inkompressiblen 
Strömung für ein vorgegebenes ebenes Schaufelgitter; VDI- 
Forschungsheft 447, 1955. 


Ein Berechnungsverfahren zum Entwurf von Schaufelgitterpro- 
filen; VDI-Forschungsheft 442, 1954. 


Verfasser: Dr. Jan Poläsek, Forschungsinstitut für 
Wärmetechnik, Praha I, Husova 8, ÖSR. 


BUCHBESPRECHUNGEN 


L. Holzer, Zahlentheorie, TeilI. VI + 2028. 
Teil II. V + 12758. Leipzig 1958/1959. B. G. Teub- 
ner Verlagsgesellschaft. Preis geb. 9,75 DM bzw. 
9,— DM. 


Diese lebendig geschriebenen und persönlich ge- 
haltenen Bändchen stellen eine Einführung dar in die 
Theorie der algebraischen Zahlkörper. Die Tendenz 
ist dabei, unter Ausnützung aller Eigenheiten dieser 
Körper möglichst schnell zu den beherrschenden 
Hauptsätzen der Theorie zu kommen, so daß allge- 
meine idealtheoretische oder bewertungstheoretische 
Gesichtspunkte nicht gebracht zu werden brauchen. 
Aufgelockert ist der ganze Stoff durch Behandlung 
vieler spezieller Typen und Zahlenbeispiele, die 
großenteils nicht nur zur Illustration der Anwendungs- 
möglichkeiten der allgemeinen Theorie dienen, sondern 
auch selbstständiges Interesse besitzen. Für den An- 
fänger störend ist eine ganze Reihe von Druckfehlern, 

und manche Sätze und Definitionen hätten etwas 
sorgfältiger formuliert sein können. Z. B. werden die 
wichtigsten Idealoperationen mit Hilfe von speziellen 
Erzeugendensystemen definiert, ohne daß auf die 
dann notwendigen Invarianzbetrachtungen hinge- 
wiesen wird. Die beiden ersten Abschnitte bringen die 
elementare Zahlentheorie mit der Theorie der linearen 
Kongruenzen, der Darstellung von Zahlen als Sum- 
men von zwei und vier Quadraten und ähnliches sowie 
einiges aus der T'heorie der höheren Kongruenzen, 
insbesondere einen Beweis des quadratischen Rezipro- 
zitätsgesetzes, bei dem Gaußsche Summen in Galois- 
feldern betrachtet werden, deren Theorie soweit, wie 
hier erforderlich, hergeleitet wird. 


Die allgemeine Theorie der algebraischen Zahl- 
körper ist dann der Inhalt des dritten und vierten Ab- 
schnittes; sie beginnt mit dem Begriff der ganzen 
algebraischen Zahl, der Ganzheitsbasis und ihrer Dis- 
kriminante, dem Fundamentalsatz der Primidealzer- 
legung, fährt mit Anwendungen des Minkowskischen 
Gitterpunktsatzes auf die Klassenzahl und auf die 
Diskriminanten fort und kommt nach der Theorie 
der Einheiten und der Körperdifferente zur Theorie 
der Relativkörper. Der allgemeine Dedekindsche 
Diskriminantensatz erscheint hier als Folge der 
Hilbertschen Theorie der relativ-Galoisschen Körper, 
die bis zur Struktur der Faktorgruppen von Zerle- 
gungs-, Trägheits- und Verzweigungsgruppen durch- 
geführt wird. 

Der letzte Abschnitt schließlich, der der Vorberei- 
tung der Klassenkörpertheorie dient, bringt zunächst 
die transzendente Bestimmung der Klassenzahl ein- 
schließlich der expliziten Formeln für quadratische 
Zahlkörper. Die ambigen Idealklassen relativ-zy- 


klischer Körper von Primzahlgrad und deren Einheiten- 
hauptgeschlecht werden mit Hilfe eines allgemeinen, 
bisher noch nicht publizierten Isomorphiesatzes von 
Rella behandelt. Den Schluß bilden Sätze über die 
Klassengruppe gewisser Kreisteilungskörper und eine 
Verschärfung für die Schranke der kleinsten Normen 
in den Klassen reeller quadratischer Zahlkörper. 


Berlin H. Reichardt 


G. Doetseh (o. Prof. a. d. Universität Freiburg) 
Einführung in Theorie und Anwendung der 
Laplace-Transformation. Ein Lehrbuch für 
Studierende der Mathematik, Physik und Ingenieur- 
wissenschaft. (Mathematische Reihe, Band 24.) 301 S. 
m. 40 Abb. Basel und Stuttgart 1958. Birkhäuser 
Verlag. Preis geb. 39,40 DM, brosch. 23,40 DM. 


Seinen bekannten umfassenden Standardwerken 
über die Theorie und die Anwendungen der Zaplace- 
Transformation fügt der Verfasser nunmehr ein ein- 
führendes Lehrbuch über denselben Gegenstand hin- 
zu, das für Studierende der Mathematik, Physik und 
Ingenieurwissenschaft gedacht ist. Die Art der Dar- 
stellung ist dem Zweck, den sie verfolgt, hervorragend 
angepaßt. In kluger Beschränkung wird von der 
Theorie nur das gebracht, was für die Anwendungen 
als Grundlage wichtig ist, das aber mit großer Klar- 
heit und mit aller wünschenswerten und doch nicht 
übertriebenen Ausführlichkeit. Wenn dabei — auch 
schon bei manchen mathematischen Definitionen — 
die Verbindung zur physikalischen Vorstellungswelt 
hergestellt wird oder sogar von dieser ausgegangen 
wird, so ist dies vor allem vom Standpunkt der In- 
genieure und Physiker nur zu begrüßen. Freilich ist 
das Buch nicht etwa ein Rezeptbuch. So mancher, 
der die Laplace-Transformation routinemäßig hand- 
habt und sie für entsprechend einfach hält, wird viel- 
leicht enttäuscht sein, daß zu ihrer einwandfreien Be- 
gründung so viel Theorie erforderlich ist. Das Buch 
wird ihn — hoffentlich — auf verhältnismäßig ange- 
nehme und lebendige Weise von der Notwendigkeit 
dieses Vorgehens überzeugen. Daß es dazu eines ge- 
wissen Umfangs an Kenntnissen aus der Integrations- 
und der Funktionentheorie bedarf, soll nicht ver- 
schwiegen werden. 


Bereits dem ersten bis zur Abbildung der Integra- 
tion, der Differentiation und der Faltung reichenden 
Paragraphen folgt, noch bevor auf die komplexe Um- 
kehrformel eingegangen wird, eine eingehende Be- 
handlung der Verwendung der Laplace-Transfor- 
mation bei der Lösung von gewöhnlichen linearen 
Differentialgleichungen und Differentialgleichungs- 


r E Be 


‘ 


Es folgt eine eingehende Erörterung der komplexen 


Umkehrformel und ihrer Handhabung (Deformation p: 


des Integrationsweges, Residuenrechnung). Besonders 


begrüßenswert ist, auch vom Standpunkt der Anwen- 


dungen, die große Aufmerksamkeit, die der Verfasser 
dem asymptotischen Verhalten der Bild- und der 
Originalfunktion widmet. Den Abschluß bilden je ein 
Paragraph über gewöhnliche Differenzengleichungen 
mit Polynomkoeffizienten, über partielle Differential- 
gleichungen (Wärmeleitungs- und Telegraphenglei- 
chung) und über Integralgleichungen 1. und 2. Art 
vom Faltungstyp. 


Wer die übrigen Werke des Verfassers kennt, wird 
ihm für die viele Mühe Dank wissen, die er darauf 
verwendet hat, sich auf den anderen Leserkreis ein- 
zustellen, und wird überrascht sein von der Fülle 
neuer und reizvoller Wege, die er dabei gegangen ist. 


Dresden H. Heinrich 


E. Martensen und K. v. Sengbusch, Numerische 
Darstellung von ebenen und rotationssym- 
metrischen transsonischen Düsenströmungrn 
mit gekrümmtem Schalldurchgang. (Mitteilun- 
gen aus dem Max-Planck-Institut für Strömungsfor- 
schung und der Aerodynamischen Versuchsanstalt, 
Nr. 19.) 112 8. m. zahlr. Tabellen und Abb. Göt- 
tingen 1958. Selbstverlag Max-Planck-Institut für 
Strömungsforschung. Preis brosch. 11,— DM. 


Zur Berechnung von Lavaldüsen braucht man als 
Anfangswerte zum Einsatz des Charakteristikenver- 
fahrens sehr genauen Aufschluß über die Geschwin- 
digkeit nach Richtung und Größe in der Umgebung 
des engsten Querschnittes. Üblicherweise arbeitet 
man hier nach der Meyerschen Reihe, wobei Glieder 
bis zur 6. Ordnung berücksichtigt werden. Bekannt- 
lich konvergiert diese Reihe nur sehr langsam. Die 
Verfasser haben daher unter Benutzung einer elektro- 
nischen Rechenmaschine die Reihenentwicklung 
weitergetrieben bis zu den Gliedern 23. Ordnung für 
das Potential und sowohl für ebene als auch rotations- 
symmetrische Störungen die Koeffizienten bestimmt. 
Das umfangreiche Tabellenmaterial enthält als Para- 
meter-die Isentropenexponenten x = 9/7; 1,3; 1,4 und 
5/3. 

Mit Hilfe dieser Werte gelang die Durchrechnung 
von Düsen so gut, daß die Kontrolle mit der Konti- 
nuitätsbedingung Fehler ergab, die kleiner als 10/go 
waren. Praktisch wird sich diese Genauigkeit sicher- 
lieh nicht erreichen lassen, da der Einfluß der Grenz- 
schicht sich dann geltend macht. Als Beispiele wer- 
den zwei ebene Düsen vorgelegt, die auf M = 1,75 
und M = 2,50 expandieren. Bei den Entwürfen wurde 
gefordert, daß die Konturen keine Krümmungssprün- 
ge aufweisen. 


Dresden W. Albrine 


E. Becker, Reibungswirkungen beim Rohr- 
windkanal. (Mitteilungen aus dem Max-Planck- 
Institut für Strömungsforschung und der Aerodyna- 
mischen Versuchsanstalt, Nr. 20.) 748. m. 24 Abb. 
Göttingen 1958. Selbstverlag Max-Planck-Institut für 
Strömungsforschung. Preis brosch. 8,— DM. 


Da bei Überschallgeschwindigkeit der Energiebe- 
darf eines stationären Windkanales hoch liegt, hat man 
frühzeitig angefangen, unter Ausnutzung des Speicher- 
prinzips Kurzzeitversuche durchzuführen. Eine sehr 
einfache Ausführungsform ist der sogenannte Rohr- 
windkanal, der aus einem zylindrischen mit Druck- 
luft gefüllten Rohr besteht. Das eine Ende bleibt ge- 
schlossen, das andere Ende mündet in eine für paral- 
leles Ausströmen ausgelegte Lavaldüse, deren Münd- 
dung vor Versuchsbeginn mit einer Membran ver- 


y n „ ar er Er IE 
systemen mit konstanten Koeffizienten und beim An- 
lese der linearen Differentialgleichung. 


Modell von 16 m Länge und 30 mm Inner 
ser nicht wesentlich größer war. Erwähnt 
sätzlich geprüft wird ein Vorschlag, di 
lichen Ruhedruckabfall infolge der Reibu 
pensieren, indem das zylindrische Rohr durch ein ] 
mit veränderlichem Querschnitt ersetzt wird. 
Sehr eingehend wird die Theorie des Vorganges I 
gelegt. Eine Impulsgleichung für instationäre k 
pressible Grenzschichten mit den i 
dieser speziellen Strömung wird integriert. 
nicht reibende Strömung wird eine nichthom 


Randbedın 


ogene 


Wellengleichung aufgestellt, wobei Wärmezufuhr zu- 
a 


gelassen ist. 


: . 
Dresden _W. Albring 


F. Smithies, Integral Equations. (Cambridge 


Tracts in Mathematics and Mathematical Physics, 
N0.49.) X + 1728. Cambridge 1958. Cambridge Uni- 
versity Press. Preis geh. 27s. 6d. 


Durch diese neuerschienene Monographie wird die 
Literatur über Integralgleichungen wesentlich berei- 
chert. Verf. gibt auf relativ kleinem Umfang eine er- 
schöpfende Darstellung im wesentlichen über die 
Theorie der linearen Fredholmschen Integralglei- 
chung 2. Art mit regulärem Kern. Der Stil ist mathe- 
matisch modern und durch geeignet gewählte Symbo- 
lik außerordentlich präzis; dennoch verlieren sich die 
Formulierungen nicht in allgemeinen Betrachtungen 


über lineare Operatoren, wie das nicht selten in Bü- 


chern über Integralgleichungen der Fall ist. Hier ste- 
hen die spezifischen Eigenschaften der Integralglei- 
chung konsequent im Vordergrund. 


Die Einheitlichkeit der Darstellung gelingt durch 
geeignete Voraussetzungen: Integration im Sinne von 
Lebesgue, Betrachtung von L®-Kernen und Benut- 
zung des Begriffes der relativ gleichmäßigen Konver- 
genz. 

Haupthemen des Inhaltes sind: Schmidtsches Ab- 
spaltungsverfahren (determinantenfreie Sätze), Fred- 
holmsche Theorie, hermitesche Kerne, singuläre 
Funktionen und singuläre Werte. Das letztgenannte 
Kapitel beinhaltet auch eine Diskussion der normalen 
Kerne und einige wichtige Sätze über Integralglei- 
chungen 1. Art. Die für die Durchführbarkeit der 
Theorie benötigten Hilfsmittel werden ausführlich 
entwickelt (Orthonormalsysteme. Satz von Riesz- 
Fischer). 

Das meisterhafte Gelingen dieses Buches beruht vor 
allem darauf, daß der Verf. selbst maßgeblicher For- 
scher auf dem Gebiete der Integralgleichungen ist, dem 
man im besonderen die Erweiterung der klassischen 
Theorie für stetige Kerne auf quadratisch integrable 
Kerne verdankt; dadurch gelangen die Betrachtungen 
erst zu einer schönen Abrundung und größeren An- 
wendbarkeit. 


Das Literaturverzeichnis gibt Anschluß an alle be- 
nutzten Quellen. Ein Sachverzeichnis und eine Liste 
für die benutzten Symbole erleichtern weiterhin die 
genußreiche Lektüre. 

Anwendungen und Beispiele liegen nicht in der Ab- 
sicht des Buches. 


Dresden 


P.H. Müller 
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Dr. W. Traupel (o. Prof. a. d. ETH Zürich), Ther- 
mische Turbomaschinen. Erster Band: Thermo- 
dynamisch-strömungstechnische Berechnung. XII + 
407 S. m. 402 Abb. u. 6 Tafeln. Berlin/Göttingen/Hei- 
delberg 1958. Springer-Verlag. Preis geb. 58,50 DM. 


Der Verfasser legt den ersten Band eines umfassen- 
den Werkes über die thermischen Strömungsmaschi- 
nen (Dampf- und Gasturbinen und Kreiselverdichter) 
vor. Seit das Standardwerk von A. Stodola im 
Jahre 1924 zum letzten Male erschienen ist, hat der 
Fachwelt ein ähnlich umfassendes Werk gefehlt. Das 
mußte umso schmerzlicher empfunden werden, als in 
‚den seitdem vergangenen drei Jahrzehnten die wissen- 
schaftliche Forschung besonders auf stömungstech- 
nischem Gebiet eine gewaltige Entwicklung erlebt hat, 
deren Ergebnisse aber nur verstreut in den verschie- 
densten Zeitschriften veröffentlicht wurden und den 
praktischen Turbomaschinenbau nur unvollkommen 
befruchtet haben. Diese Diskrepanz zwischen wissen- 
schaftlicher Erkenntnis und ihrer praktischen Ver- 
wirklichung sucht der Verfasser durch sein zweibän- 
diges Werk zu beseitigen. 


Der vorliegende erste Band ist den thermodyna- 
mischen und strömungstechnischen Grundlagen ge- 
widmet, während dem nachfolgenden zweiten Band 
die Behandlung der Regelung, Festigkeitsfragen und 
Probleme der Dynamik, also wahrscheinlich Wellen- 
und Schaufelschwingungen vorbehalten bleiben soll. 
In dem vorliegenden Band wird nach eingehender 
Behandlung der thermodynamischen und strömungs- 
technischen Grundlagen für Gase und Wasserdampf 
unter Anwendung der modernen Methoden der 
Strömungslehre auf das Einzelprofil und das Schau- 
felgitter in ein- und mehrdimensionaler Betrachtung 
die Auslegung der Einzelstufe und der mehrstufigen 
Turbomaschinen anhand von kennzeichnenden Be- 
rechnungsbeispielen aufgezeigt. Bemerkenswert ist 
dabei, daß der Verfasser versucht, anstelle der bisher 
durchweg üblichen durch Verlustbeiwerte der einzel- 
nen Stufenelemente gekennzeichneten Berechnungs- 
verfahrens einen Radwirkungsgrad einzuführen, über 
dessen Zweckmäßigkeit der praktische Gebrauch wird 
entscheiden müssen. Das Buch wird mit der Behand- 
lung der Dichtungsfragen und der bei Anderung des 
Betriebszustandes auftretenden Probleme beschlossen. 


Das Buch ist vor allem durch drei Charakteristiken 


‘ausgezeichnet: der Verfasser versucht in allen Teilen 


miö großem Geschick das Gemeinsame der Turbo- 
maschinen herauszustellen und so zu allgemeinen Aus- 
sagen zu kommen; als zweites ist die Verknüpfung der 
theoretischen Erkenntnisse mit den konstruktiven 
und betrieblichen Fragen der praktischen Wirklich- 
keit zu nennen, und als drittes die den Berechner von 
Turbomaschinen besonders berührende Darstellung 
moderner Berechnungsmethoden als praktische An- 
wendung der entwickelten Theorien, die gleichzeitig 
die Problemstellungen und ihre Lösungsmöglichkeiten 
verdeutlichen. 

Das Buch bildet mit der Klarheit, Folgerichtigkeit 
und Gründlichkeit seiner Darstellung ein vorzügliches 
Lehrmittel für den praktischen Ingenieur wie für den 
Studierenden, setzt allerdings beträchtliche Vorkennt- 
nisse auf mathematischem, hydro- und aerodynami- 
schem Gebiet voraus. Dem Mathematiker und Me- 
chaniker selbst bietet es eine Fülle interessanter An- 
regungen für die praktische Anwendung physikali- 
scher Gesetzmäßigkeiten auf die Vorgänge in den 
thermischen Strömungsmaschinen. 


Dresden H. Faltin 


Dr. P. Lorenzen (o. Prof. a. d. Universität Kiel), 
Formale Logik. (Sig. Göschen, Bd. 1176/1176a.) 
165 S. Berlin 1958. Walter de Gruyter & Co. Preis 
brosch. 4,80 DM. ; 

In diesem Bändchen wird der Versuch unternom- 


men, eine konzentrierte Einführung in das Gesamt- 
gebiet der formalen (mathematischen) Logik zu geben 


Ein erstes Kapitel behandelt die aristotelische Syllo- 
gistik. Es folgen drei (!) Kapitel über Aussagenlogik. 
In einem fünften Kapitel wird der Prädikatenkalkül 
der ersten Stufe ohne Identität und in einem sechsten 
und letzten Kapitel wird der Prädikatenkalkül der 
ersten Stufe mit Identität behandelt. Die Prädikaten- 
kalküle höherer Stufe werden ‚als nicht zur formalen 
Logik gehörig‘ aus der Logik verbannt. Die vor- 
liegende Darstellung zeichnet sich dadurch aus, daß 
fast alle auch nur einigermaßen üblichen Bezeich- 
nungen und ‚Symbole durch neue ersetzt werden. 
Noch schlimmer ist aber, daß auch sachlich in jeder 
Hinsicht von der bisherigen Entwicklung abgewichen 
wird, so daß selbst der Fachmann auf dem Gebiet 
der mathematischen Logik viele Dinge, die zum all- 
gemeinen Bestand der mathematischen Logik ge- 
hören, nur mit großer Mühe wiedererkennt. Das geht 
so weit, daß viele grundlegende Entwicklungslinien 
vollständig ignoriert werden. So wird z. B. in dem 
Kapitel über Syllogistik die heute wohl allgemein an- 
erkannte Interpretation der aristotelischen Logik 
durch Lukasiewiez mit keinem Wort erwähnt; den 
umfangreichen Ausführungen zur Aussagenlogik fehlt 
der geringste Zusammenhang mit den fundamentalen 
Untersuchungen der Warschauer Schule, usw. (In- 
teressant ist in diesem Zusammenhang, daß z. B. der 
Enzyklopädie-Artikel über mathematische Logik von 
Hermes und Scholz im Literaturverzeichnis nicht 
zitiertist; diegrundlegenden Lehrbücher vonHilbert- 
Ackermann, Hilbert-Bernays und Church sind 
zwar im Literaturverzeichnis- aufgeführt, es wird im 
Buch selbst jedoch an keiner Stellean die dortigen Ent- 
wicklungen angeknüpft). Daher kann die vorliegende 
Darstellung zwar dem Kenner einige neue Aspekte 
vermitteln, ist aber leider kaum dazu geeignet, dem 
Leser, der die moderne Logik erst erlernen möchte, 
eine auch nur einigermaßen richtige Vorstellung von 
ihr zu vermitteln. 


Berlin G. Asser 


W.V.D.Hodge, Theorie und Anwendungen 
harmonischer Integrale. VIII + 247S. Leipzig 
1958. B. G. Teubner Verlagsgesellschaft. Preis geb. 
15,75 DM. 


Verf.-hat die Theorie der harmonischen Integrale in 
der Absicht begründet, eines der fundamentalen Er- 
gebnisse der klassischen Mathematik, den Riemann- 
schen Satz von der Existenz überall endlicher Inte- 
grale auf einer Riemannschen Fläche, in allgemeinere 
Zusammenhänge einzuordnen. Die 1941 erschienene 
erste Auflage des Originals enthält eine Zusammen- 
fassung der bis dahin auf diesem Gebiet erschienenen 
Arbeiten mit dem zentralen Satz von Hodge über die 
Existenz harmonischer Integrale mit vorgeschriebenen 
Perioden auf einer Riemannschen Mannigfaltigkeit. 
Die 1952 erschienene zweite Auflage unterschied sich 
von der ersten nur in wenigen, aber wichtigen Korrek- 
turen und Abänderungen; die sonst seitdem erzielten 
weiteren Fortschritte wurden jedoch nur in ein Li- 
teraturverzeichnis aufgenommen. Diese Fassung — 
ohne weitere Änderungen oder Hinzufügungen — 
liegt der Übersetzung zugrunde, die für die weitere 
Verbreitung dieser tiefliegenden Theorie wesentlich 
beitragen wird, 

Der ganze Stoff wird in einer sehr gut lesbaren 
Form dargeboten. Die beiden ersten Kapitel stellen 
die später benötigten Vorkenntnisse zusammen, und 
zwar enthält das erste die lokale und die globale The- 
orie der affin zusammenhängenden und der Riemann- 
schen Mannigfaltigkeiten, also deren Differential- 
geometrie und Topologie, während das zweite den 
alternierenden Differentialformen und ihren Inte- 
gralen gewidmet ist und in den beiden Sätzen von 
deRham gipfelt, deren erster die Existenz von p- 
Formen mit gegebenen Perioden sichert und deren 
zweiter besagt, daß die geschlossenen Formen (For 


_ men mit verschwindendem alternierendem Differen- 


tial), deren Perioden null sind, Nullformen, d.h. 


‚alternierende Differentiale geeigneter Formen sind. 
Im dritten Kapitel wird eine alternierende Differen- 


tialform als harmonisch auf einer Riemannschen Man- 

nigfaltigkeit definiert, wenn sie zugleich mit ihrer du- 

alen Form geschlossen ist. (Im klassischen Fallkommt 
man so auf die Cauchy-Riemannschen Differential- 

gleichungen.) ; 

Der von Hodge aufgestellte Hauptsatz der ganzen 
Theorie besagt nun die Existenz harmonischer Inte- 
grale mit gegebenen Perioden. Zu seinem Beweis 
wird die Fredholmsche Theorie in ähnlicher, jedoch 
wesentlich komplizierterer Weise wie in der Poten- 
tialtheorie herangezogen. Die beiden letzten Kapitel 
enthalten Anwendungen, und zwar beschäftigt sich 
das vierte mit algebraischen Mannigfaltigkeiten. Da 
hier Invarianten nur gegenüber birationalen Trans- 
formationen gesucht werden, ist es auf viele Weise 
möglich, eine solche Mannigfaltigkeit mit einer Rie- 
mannschen Metrik zu versehen. Wendet man daher 
die Ergebnisse des dritten Kapitels auf eine solche 
Mannigfaltigkeit an, so muß man sich nachträglich 
wieder von der Metrisierung zu befreien suchen, und 
es erweist sich dann eine gewisse Periodenmatrix 
harmonischer Integrale als birational invariant. Das 
fünfte Kapitel schließlich bringt Anwendungen auf 
die Theorie der kontinuierlichen Gruppen. Auf einer 
Mannigfaltigkeit wird durch eine halbeinfache, ab- 
geschlossene und transitive Transformationsgruppe 
eine Metrik erzeugt. Die zugehörigen harmonischen 
Integrale liefern dann einen besonders durchsichtigen 
Zugang zu gewissen von Cartän untersuchten Grup- 
peninvarianten. 


Berlin H. Reichardt 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


mathematik der GAMM. Entsprechend der spe 


_ Begrenzung betrachtet. Auch der genauer 


8. 11. 1957, zusammengestellt v 
m. 121 Abb. München 1958. R. 


Preis geb. 16,80 DM. Bar er 
Den Einsatz moderner Rechenanlagen in de 
lungstechnik zu fördern ist das Hauptziel 
Buches. Es enthält Vorträge und Diskussions 
von einer Fachtagung des Ausschusses für R 


Struktur der Aufgaben der Regelungstechnik s 
Verwendung von Analogiegeräten im Vordergrund. 

Wie die manchmal nicht ganz einfache Wahl der 
Skalenfaktoren vorzunehmen ist, wird in einem Vor- 
trag von Herschel gezeigt. Analogiegeräte sind auch 
besonders geeignet zur Untersuchung des Einflusses 
vorkommender Nichtlinearitäten. Ein V. von 
Bühler behandelt z.B. Regelsysteme mit Reibung, in 
einem Diskussionsbeitrag von Ernst wird die Frage 
der Wirkung von Nichtlinearitäten wie Totzeit oder 


allgemeinen langsamer arbeitende Digitalrechner hat 
Bedeutung in der Regelungstechnik, er eignet sich 
z.B. gut zur Untersuchung des Stabilitätsverhaltens 
(Vortrag von Bukovics). 

Eine ganze Anzahl weiterer Vorträge gibt Einblick 
in einige in der Regelungstechnik auftretende Pro- 
bleme und die Art ihrer maschinellen Behandlung. 
Die Art der Zusammenstellung aus Einzelvorträgen 
gestaltet das Buch sehr lebendig und wird ihm sicher 
einen weiten Leserkreis gewinnen. 


Berlin K.-H. Bachmann 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


M. Schuler, Mechanische Schwingungslehre, 
Teil II: Mehrfache Schwinger. 1./2. Aufl. V-+ 
150 S. m. 100 Abb. Leipzig 1959. Akademische Ver- 
lagsgesellschaft Geest & Portig. Preis geb. 14,— DM. 


Sammelband zu Ehren des 250. Geburts- 
tages von Leonhard Euler. Verantwortl. Redak- 
tion K. Schröder, Dtsch. Akademie d. Wissenschaften 
zu Berlin. X + 336 S. m. 20 Abb. Berlin 1959. Akade- 
mie-Verlag Berlin. Preis geb. 58,— DM. 


G. Matz/P. Gayer, Berechnung von gußeiser- 
nen emaillierten Druckbehältern. Messungen 
an einem 6 cbm-Rührwerksbehälter. XVI -- 233 8. 
m. 97 Abb. Berlin/Göttingen/Heidelberg 1959. Sprin- 
ger-Verlag. Preis geb. 37,50 DM. 


(. Berge, Espaces Topologiques— Fonctions 
Multivoques. (Collection Universitaire de Mathe- 
matiques) XI-- 2728. m. 47 Abb. Paris 1959. 
Dunod Editeur. Preis geb. 3.400 F. 


Abhängigkeit der Eigenschaften Hochpoly- 
merer von der Vorgeschichte des Materials. 
Vorträge und Diskussionen der 3. Tagung, Bad Nau- 
heim 22.—23.4.58. 1428. m. 149 Abb. Darmstadt 
1959. Steinkopff-Verlag. Preis brosch. 20,— DM. 


H.F.Schwenkhagen, Allgemeine Wechselstrom- j 
lehre, Zweiter Bd.: Vierpole, Leitungen, Wellen. 

XII + 4418. m. 335 Abb. Berlin/Göttingen/Heidel- 4 

berg 1959. Springer-Verlag. Preis geb. 39,—DM. ’ 

H.Schilt, Elektrizitätslehre. 216 S. m. 187 Abb. er 

f 


Basel/Stuttgart 1959. Birkhäuser-Verlag. Preis geb. 
24,— DM. 


G. B. Thomas, Elements of Caleulus and Ana- > 
lytic Geometry. X + 580S. m. 268 Abb. London 
1959. Addison-Wesley Publishing Company. Preis 
geb. $ 7.50 bzw. 57 s. / 


D.K. Cheng, Analysis of Linear Systems. 
X111 -+ 431 S. m. 272 Abb. London 1959. Addison- 
Wesley Publishing Company. Preis geb. $ 8.50 bzw. 
64 s. br 


W. Macke, Quanten. Ein Lehrbuch der theore- 
tischen Physik. XII + 4948. m. 127 Abb. Leipzig 
1959. Akademische Verlagsgesellschaft Geest & Por- 
tig. Preis geb. 29,50 DM. 


H. Kühl, Probleme des Kreuzstrom- Wärme- 2 
austauschers. VII+83S. m. 34 Abb. Berlin/ a 
Göttingen/Heidelberg 1959. Springer-Verlag. Preis 
geb. 10,50 DM. 
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